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CHAPTER  I 
CHXPTIO0  0?  QDLID  SELIUM 


flare  Cas  Colldo 

physical  properties  of  the  solid  phases  of  'the  'rare  gft»H 
elecients  are  of  particular  Interest  because  their  atomic  electrons  occupy 
fully  closed  sheila,  and  the  only  chaBical  binding  force  is  the  weak  van 
der  t&als  interaction.     3he  physical  properties  of  these  eolids  are  thus 
in  aany  respects  sore  open  to  analysis  than  nost  other  solids.    Further, 
the  veafcaeas  of  the  iaterpsxticle  force  gives  rise  to  .te  quo. 

©echenioal  effects  in  the  aotion  of  the  atono  in  the  censed  state. 

Sobbs  and  Jones    have  written  an  extensive  review  of  the  e^peri- 
nanta!  and  theoretical  work  on  these  solids,  with  enphaais  c  >?,rgon. 

2&ey  show  the  importance  of  the  q^ituarfieciianical  aero-point  energy  by 
listing  the  ratio  of  aero-point  energy  to  aches  ojqgy  for  the  nn 

^aseo  .yr&ev  of  decreasing  ataalc  weight);     MM*  0.031; 

krypton,  Q,QJk$  arigoa,  0.0975  neon,  0.311;  b*Hu 

ordinary  classical  lattice  dynamics  applies  to  xenon,  krypton,  and  arson, 
with  small  qjuttituaHoechanicel  corrections.      ban  can  be  treated  adequately 
by  including  aahawaonie  tanas  in  the  lateratoaic  potential.    Eowev<- 
heliua  is  of  such  low  atoraic  $ass  Urn    I     o  powerful  techniques  in 


*!•  IU  Sobbs  and  C.     .    .'ones,  2ieory  and  Properties  of 

AHjHfe     .seporto  on  Progress  in  Hr/oicn,  Vol.  XX,  5i6  (1957)* 


be  BB&bspeA*      Further,  there  is  also  an  isotope  of  atoxaic  wight  three* 
vshlch  has  spin  IJ2  sod  obeys  the  Feral  -S&rnc  static 

IMs  thaaia  is  directed  to  the  task  of  d&v    ■  ,he  grouad  state 

^  la  a  vny  <wbioh  iaolu&ee  qnanttsa  effects  la  Use  basil 
'*he  description  of  the  abound  state  «11X  be  need  in  tuv  .alculr 

physical  ojtjaatitiea  *&ioh  will  be  coopered  with  e;g»ivisi©ntaI2y  observed 


Bernardeo  and  Priaakoff3  (hereafter  termed  B.P.)  have  prose- 
's 
a  theory  of  solid  BeJ  end  have  drawn  extensive  conclusions  tan  their  ep- 

proxiiaate  grouad»ata&e  wave  X'unctlon.    She  present  wosfc  is  ssali:. 

attest  to  give  a  raore  accurate  description  of  the  ground  than  the 

ftoanlatlcxi       §*?< 

Properties  of  Atoalc  and  Solid  Belim 

II    »      II     iWllll     nil I        II  III.  ■!...«  in       .id, Hi     »    1   I   11,   ,|,   ,l,.lll»    ,.«.,       1,1,11 .,».  1  I.    If  — 

k 
Bolid  He  possesses  no  such  starUiag  properties  as  the  liquid 

phase.  2here  has  been  no  such  intensive  application  of  pove.       .t?t- 

ical  tools  to  the  solid,  as  there  has  been  to  the  liquid.  However,  n 

experimental  voi&  has  been  performed,  end  theoretical  explanation  isada  of 

k 

the  observed  reoults.  Several  reviews  of  the  properties  of  solid  Zfe 

ejk,a 

have  been  published.   '  " 


a 
C.  Itoub  and  J.  B«  DMBdolej    Solid  lieliun.     Procreoo  in  lov 
fae^rature  fflyelcs,  Vol.  II  tfj.  J.  Gortcr,  ed!   •        Xatarseleac"  .      s, 

Hew  Yorli ,  M?5  p^  3: 

%U  aernardes  and  &  lYtmimnfr      Ehy.         . .  3J^>  96©  (1961). 

;.oaaoa,    Superflulds,  Vol.  21  ( Wiley  and  Sons,  Sev  lor  1%)* 

fcti&iaa,     liquid  Beliuci  (Coabxlclee  University  Press, 
'.bridge*  1959  )• 


-3*" 


$,  Eer  has  beer,  aside  available  in  increasing  eoour.. 
end  there  has  Ixsen  considerable  sracent  u  th  this  isotope.    3road 

discussion  or  the  properties  of  liquid  and  solid  lie*'  ic  contained  ia 

Toe  nuclear  aoaent  of  He*'  is  -  8*  lear  an©-  Its 

atoaic  nemo  (on  the  0      scale)  is  3-01^3*  &&1*  that  of  lis    is  U.OO 


Properties  uhich  are  directly  due  to  the  atonic  electrons  are  the 
for  the  two  isotopes  aside  frees  an  oi-rfcrosaeiy  mall  *ive  saass  oorae 

tion.     -23©ae  effects  include  atomic  spectra  and  the  interatoaic  p..  i» 

At  ataospberlc  pressure,  He";    i        io  ,  :.xjraaina  liquid  at  absolute 
aero  imperative*     Crystallisation  teho©  place  on  application  o 
30  eta  af  pressure.    A  pressure  of  tvem:.  .8  needed  to  solidify 

He1*. 

1 
Ihe  properties  of  solid  &r  havo  not  been  thorougily  explored 

earperiaentally.    Much  nuclear  resonance  er-pcriaaatation  has  bee  cd 


■- 

jB.  F.  Ilbaael,    3ao  lov  ^eaperatvi  -.wtlao  of  Heliv 

Progress  in  lay  Tenpemture  ]Physics»  Vol.  I  (C.  J.  Gorte?,  editor,  lntts?~ 
science  Publisher,  stew  tforfc,  19o*?)#P»  T 

'V.  P.  fteshUov  and  K«  H.  Sinoveva,     ifcpcrts  on  Progress  in  Physics, 
Vol.  BIX,  90*  (1#59). 

(Chio  State  Uhiversit^,  IS  57). 

-oliua  flareo  (J*  0.  Daunt,  editor,  Chio  State  University  ?wm»o, 
Columbus ,  Ohio,  ii'o35'. 

i 

proceedings  o  .  ... ,.  .■•  ^tionol  Conferen . 

Bsmti-    ..v?  Physics  [0.  K,  arenas  ana    •     .    bills  hUirt)  editors, 
flsiversity  of  55oronto  l^ress,  1961),  p.  563. 

13.  Ti.  Cxllfy  and  £.  F.  ltaaac-  ,  '  (I  and  :io  .     Progress 

■a  Urn  Smpe        -.-  Hayaics,  Vol.  XX3      -    '■•  Porltrj  c.i-U"..-,  to 
liahed}. 


~«v 


out,    W*  H.  ^rbani:  has  i*eviewed  this  «  '        Data  on  the 

relationships  near  the  aelti&g  curve  haw  been  grfolish*  .iy> 

IS 

0,  0.  clydoriak,  end  &  X*.  IHUe/  and  by  Grilly  and  Mills. 

theraal  condmtivity  have  been  ptioHefced  by  !•  J.  t&Oker  and 

Ochueh,  Crilly  and  MiBs,  £h  X-ray  e^fftactioa,  have  found 

that  solid  lie    exists  in  two  phases.     9m  scrolled    a    phase,  a  bo&y- 
center  cubic  configuration,  CJ&sta  at  pressures  below  about  300  ate*  and 
-Use    0    phase,  of  hexa&aoftl  cloee-paaicad  ■IrmlWHj  exiate  at  lugh' 
pressures.     Ihis  phase  transition  was  also  observed  by  Griily  and  t&Us* 
as  a  systematic  discontinuity  in  the  specific  voluae. 

Pujliiiinary  EjeaflurcBoate  af  specific  heat  of  the  solid  have  bean 
nade  by  Bayards,  UauEi,  Brewer,  flaunt  and  tfcUilliasr-  The  pressure* 

tespera&urc  dependence  of  the  taaltine  ourve  Jias  been  obtained  to  0« 
by  ti»  smae  group.     2hey  find  the  entropy  to  have  o  constant  value  close 
to    ft&&    over  tlio  range  0,0? 

She  principal  anooaly  in  the  behavior  of  the  solid  lies  in  a  wiui- 
EJUEi  in  the  Belting  pressure  curve,  plotted  as  a  function  of  tee:  re* 

at    p  •  8£»3  eta*  *  »  »Jft*  ft.    She  a  oa  of  the  Btnirara  is  aeaoapen:< 

by  interesting  ;  lading  a  aa@atiw  thenaal  w-ganai  mt, 

and  the  consetr  ..saibility  of  oeltine  the  solid  by  lowering  Itf 

temperature  at  constant  x:aft?soure.     r  'iapeyrca  equate'. 


13 
TO  (1. 


L23S.  n.  Gv  ..nd  IU  L  IHllo,     Aon.  Phys.  Q,  1,  (2&S9) » 

..:&,  s.         .illy,  am.:  Us,    Bye.  fern  up, 


""lie. 


2u*  e:rlet*oo©  of  higher  thcraodyaeiaic  ovte?  la  the  liquid  than  in  the 
solid  vouid  yield  ouch  behavior 

She  neliua  Interpartiele  Potential  gad. 
'fee  Many  a*fr  Hoailtanian 

Ihc  stationary  states  of  an  aggregate  of  bail -in  atoas  are  solution* 
of  the  Gchrodingar  equation 


-  ,. :_  ,  (i) 


R^  f  „  ■ 


vhere  tl.      Is  the  nany  body  HeE&itonlaa  involving  all  the  nuclei:!  and  all 

the  electron©  of  the  solid* 


S 


«  •  £  &-      .2!—       ~  £'  2 

i*l^e  *   a-1®*  'a     1*4  riJ 

*  Val(Xi  .  .  .  ^  yx  .   .  .  %)  •  Vm{yi   .  ,  .  y,).     (2) 

Here  1  labels  one  of  the  K  elac'     a  labels  one  of  the 
8  nuclei*.,  the  K±     label  all  the  coordinates  (includinc  epln)  of 
electron  i.   tfuelaor  aoordinates  are  described  by  the  yfs»   Vr.  .  tmA 
V   are  the  nueleiua^iaetron  and  mjcleua~nueieus  potential  energy. 

lowing  the  adiabatic  approximate  of  Sora  and  Oppenhein    r  wi 
separate  if    into  a  port  depending  upon  electron  coordinates  that  contains 


3ora  and  J.  &•  Opsponheioer,  Ann.  Hjyai!*,  2b,   fcg     67). 


nuclear  coordinates  go  pasacieters,  end  &  pert  aapasdine  only  on  nuclear 
-.dinatea 

i       •  •  •  ^  fv ■  ■  -  V  0  (*i>  ■  *  *  %> 

We  dei'ine  the  part  of    IJ^    defending  on  ate  coordinates 

end  nuclear  chorees  as    H  .     Ihea 

(3) 
vhere  we  ohov  the  f>chr6din©er  equation  for  the  electrons,  from  "tihiob  v<*? 

obtain  the  electron  energy  and  state,  for  a  particular  arrangement  of  the 

nuclei!  described  by  the  valuee  of  the  parameters    (j.j  y    .   .   •  ;? 

15 
through  application  of  the  adiabatic  appraadjaatlon,  '  -where  ve 

lali  the  remillilll  dopeadencc  of    Ee    end    $    on  the    y'e,    aestiaing 


only  tJiat  the  electron  motion  can  adju  notant  to  the 

positions  of  the  nuclei!,  and  vlth  the  notion  atiil  deierolned  by  (3), 
vo  obtain 

z  - §  V*  1  +  Ee(yr  H  .  . ,  rjt  -  #  CM 

a*i 

We  see  fron  thie  dovolopaent  that  the  in'..    iear  p 
{for  that  is  the  role  it  playe  in  (4))  1©  a  function  of  the  pooitions  of 
all  the  ai&lglly  and  results  fron  the  presence  of  the  atomic  electrons 
and  nuclear  charges. 


lei 

yA  coag?iete  development  of  the  adiabatic  appro^lnatioa  ■ 
la  the  prase:.  >  by  P.  Seita  in  M.  heory  ..Ids 

(ncQrmt  Bill  Book  Co.,  Inc.,  Ifev  Torte,  19^0)  p/W.' 
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*2a  obtain  an  explicit  torn  detailed  study  of  (k)  is 

required. 

mrsdafelder,  Curtiao  and  Bird     have  described  the  caaraacOy  used 

phenciaenoiogicai  forne  o:°  the  Interne  :.el  of  M    atcras,  and 

17 
have  given  an  overview  of  the  theory.  .      has  soispared  several  of 

18 
these  foitae  of  the  potential.    Moore      has  nade  a  recent  0  butiao 

the  theory  of  the  repulsive  core  of  the  potential. 

At  distances  large  coopered  to  the  Bohr  radius,  the  into 
of  two  belivia  ntoos  as  described  by  (k)  is  characterised  by  the  van  cu 
Waals    .  -ceractlon*  *&ieh  arises  froci  ftttrs  between  the  e&e 

dlpola  eoaents  nutually  induced  in  the  electron  charge  die 
Shis  interaction  energy  depends  inversely  on  the  sixth  power  of  the  oep* 
aration  distance    r    of  the  atone.    99m         1  also  higher  order  teaasj 
the  next  la  the  dipole-quedrupoie  energy*  uhlcb  is  proportional  to    r    . 
At  close  di  i  of  approach,  of  tfee  or  two  Bohr  radii,  the 

potential  is  steeply  repulsive,  roughly  ional  to    0  in  A). 

Sue  electron  clouds  are  stabilised  9  each  nucleus  ■  .^sigJy  repel 

each  other  vhon  they  start  to  overlap  because  of  both  Couksab  repulsion 
and  the  Pauli  principle* 

•She  phononeaologioal  foists  of  the  potential  have  been  obtained  by 
fitting  the.  theoretical  forsaa  for  Hie  potential  e  tengjersture  depenO- 


J.  0*  liirschfelder,  C.  F.  Curtice  WA  •  Bird,    Ifolecular 

ay  of  Oases  (Wiley  and  Bout*  Zlev  Ja)  pp."  "2&6k  and 

17J.  DeBoer,     Physics  &  ). 

NMMftj    J.  Ck«b  Phys.  ]},  **?!,  {i960). 


ence  of  the  second  virial  c  -.'it  of  the  ra#e  gcneo  In  the  gaseous 

etate,    A  ft-  the  potent!?,  b  hen  bean  calculated  with  low  scape 

ature  data  end  which  has  been  found  t  the  vlrlal  laient  v; 

.pod  .««*,  to  «  *«.;  **-.  «««  1  MM, 

hereafter  tamed  &.22.M.    !2aiij  potential  hats  the  torn: 

ft*)  -  [r?0  e"  ]  erg         r  «  a.6l 

r ' 

(5) 

V(r)  -  [977  *  -  ^        "  •  ]  **S  -61  A, 

r 

Analytical  appllen  -iaia  fom  would  "be  difficult.     However,  the  two- 

paraaeter  potential  of  lanaerd  -Jones 

9C9>  **«[Cf»^-  if)6]  <6> 

with    c  »  30-  .,    *  *  8.      •  A,    Gives  ft  very  close  1  the  B.SI. 

potential*    It©  behavior  in  the  region  of  severe  refusion  varies 
that  or  toe  :  ae  fl.K.U.  ft  at  it  will  be  apparent  later  that 

>g  deviation  takes  place  In  a  region  ft  high  repulsive  potential 

aa  to  leave  physical  i,«esulta  unofffceta.  . 

"Zhree  helita  atosao  in  prox&Bity  vat  ■  three-wey  Induced  dipole 

interaction,  i&ich  follows  the  ralatloodhi;  ■ 


•^JU  A.  Buc;a-\$icn,  J.  Hem:.  aid  II-  •  fttoaaey,     f*  EQT* 

.  (London)  A17,  .). 

aa 

2t  ia  later  den  ":*>d  that  the 

repulsive  ia  to  reduce  the  probi  penalty  obtained  fraa  the  M 

body  wave  function  prr  .:.thar  torn  of  12ie  potential 

cijuaiiy  effective  la  this  saga; 


2I?.  ifcaen,    J.  Chen.  Hsya*  21,  1007,  (1 


-   - 


"\iO:: 


*    9.8  «",tt*pab  •  *te) , 

sab  *  %g  *    c& 

lAserc    rRb    is  the  distance  bet\**n  atma    a    cad   *    neaeiired  la  Ifchr 
radii,    E^    Is  the  induoad  triple  •  eaessy*  and    E^.     is  the  two- 

body  energy,     la  cases  of  interest  to  us  .  t2bm 

Bebc^Eab  +  ^ac  *  ^be^    i0  of  the  or<2j6:''  ^  •  ^  eo  caa  ^° 

Sence  us  are  correct  ia  stating  that    EpCy^*  •   *  »  y»)    I 

scribed  adequate^  fbr  our  purposes  by  (6),  and  the  asay-body  Ha&. 

vrittc  —lie  all  Coulorib  and  c  ©ffecta  Is 

e  \ 


art.      **■      *  a,b»l 


I  l 


2his  flsmlltonioa  is  valid  under  the  following  liberal  conditions: 
1«    'itefflperatuitj  lees  than  about  i?0* 
i     aenoity  less  than  that  >spca<lltt0  u>  a 

neighbor  distance  of  2,5  A. 


; 


ZerQ-?a&at 


JSL-L 


Ail  GeUua 


We  shell  use  the  aany-body  Houiltonian  containing  the  nuclear 
position  coordinates 


If  ■      f    .  SL      .7  ^   *  i  ■  •      \ 


(i) 


with  the  potential 


V(rlJ)  *  *  e 


*- 


■1J 


\ 


'<  4 


6. 


(2) 


between  pairs  of  helluB 

Our  purpose  is  to  solve  the  Bohri>din©er  equation  for  the  ground 
state  of  the  solid,  tyQ    and  to  obtain  the  ground  state  energy  2L  , 
given  by 


E  ^,   £,  .  .  .  vn)  .  20  ^,  % 


^ 


(3) 


where    X1    describee  all  the  relevant  coordlnateo  of  particle  1,  including 

spin,  of  any.      ty0    oust  have  the  syonetry  required  by  the  statist:;. 

3 
For  lie  ,  which  has  spin  1/2  and  obeys  FeraUZftrac  statistics-,  we  r*qp 

that 

Pe0-(-  Dpt0 


-  10- 


MQPftftt 


■ 


ufcere    P    ioea  operator  which  pexoutee  the  Indices  of  the  coordinates, 
and    p    le  the  parity  Of  this  penautation.     ~be  state  in  (4)  is  necessarily 
antiayaaetric  under  interchange  of  all  the  coordinated  of  two  particles. 
lie    obeys  2ose  statistics  end  has  a  wave  fsactlca  sycsaetric  under  ; 
chen&e  of  the  coordinates  of  "cvo  particles 

We  can  use  the  ©round  otate  to  obtain  properties  of  the  solid  at 
finite  teoperature  by  superlcgposittc  i*onon&,  and  spin  waves,  o*  applying 
the  saolecular  field  theory  of  parsjaagac 

Ihe  potential  (2)  la  characterised  by  a  weak  attraction  at 


*u  ,    a  vl***  mm  1.- m,  «d  .  etr^e  repulse*,  insiae    ,-..». 

two  paraoeters    o    and    e    can  be  asai^sed  values  to  give  good  representa- 
tions of  the  potentials  between  rare  gso  atons,  end  fair  opprosdjaatlono 
other  oysteciB,  cr^ch  aa  nucieon-nuuieor,,  aitro^en-nitrogen,  etc.    Denee  we 

can  consider  a  broad  class  of  problems  by  reducing  (1),  (2)  and  (3)  to 

t 
dineasionleoo  fora.    Following  the  discussion  of  Demardes,    we  define 


J.  Daunt 


,lH.  Bernardes,  Sbsory  of  "olid  23s3,  Belima  Ihree  (Sdited 
,  Chio  ntote  University  °reoc,  Coiuribua,  Ctolo,  VjSO)  p.  115. 


m&tiZ** 


She  iaat  ralationabip  is  &L?m  by    (l/e2)  d2/^2  [vC^o)]  ■  !>"  feL^fe)* 

-wtsere  the  l&placiaa  ie  aohexaatically  repreaanted  by  the  second  partial  de- 
rivative* 

"deine  the  definitiona  (5)#  (3)  becoewo 

c  -        *  (xx .  .  .  ^       t  i  vu^  *  •  •  v  *n  *  *  * 


v'.-JE-j. 


Bgpation  (6)  applies  to  any  pbyr..  ;/st€S3  described  by  the  two 

paranatal*  potential  (2).    She  jfayeical  properties  Of  a  particular  oysten 
are  all  contained  Iff  the  paraoater    X    and  the  fllflnflfrttm  of   &*.    Values 
of    X    3*or  various  aubator^cco  are      jseaon,  7  *  10  ^;  neon,  7  8  10    J 
He4,  Go^7;  He3,  O.yXi;  neutxvn-proton  interaction  fc  ^/(\^))f 

(rQ  *  nuclear  radius  la  fenaio,  VQ  «  potential  depth  in  ISev,  X  *  0.3)* 

She  ideal  aituatioa  would  pe  oonoral  forr,  of    ijf,  X 

aa  a  contlnuoua  pareeaeter,  but  no  such  siopli^led  solution  boa  \3oa:i  ob- 
tained for  (6). 

In  the  JiHttine  caae  i»  \Aiich    X"  <  <  1>    ^iere  the  <maatun 
aapecte  of    y    can  be  neglected,  w  obtc 


*     m 


1  5?/?? 


iu 


In  tills  case,  Hm  earrsy  o:T  the  ground  ctate  ie  the  potential 
♦aergy  lattice  aun.  She  expression  (?)  is  valid  fbr  raost  solids,  and 


- 


~13~ 


ordinary  lattice  dynaadcs,  inolndinc  quantised  sound  fields,  start  vlth 
(7)  and  Sswton's  second  lea?. 

f  larger    X.,    the  kinetic  energy  torn  of  (6)  (the  soro-polnt 
enercy)  is  on  important  pert  of  the  energy  operator.    She  theories  atf 

liquid  heli'jci  and  nuclear  natter  have  the  eoacjan  central  problea  or  the 
handling  of  the  :.irjctie  energy  tern  slnulteneously  -with  the  tvo-partieie 
potential*    Solid  Ikr  and  Ife  *"  share  this  difficulty. 

gbeorleo  of  Solid  jte^pgfag 
reveroi  author©  have  ceaaiderod  the  dynerrlae  of  solid.  Ife  . 
brief  description  of  their  .vori;  pertinent  to  our  problsEt  io  giv  t5, 

(a)    Hooton"  hoc  considered  the  ahope  of  the  potential  around  the 
lattice  site.    She  zero-point  notion  orqpande  the  solid  so  that  the  nearest 
interparticle  distance  lies  outside  the  p**«*riT>  of  the  two-particle  poten- 
tial.   Ihe  resultant  I         •  sua  shows  I  anhoraoji.  Hooton  < 
procdrmtes  the  resulting  well  tsy  a  Imrsaantc  oscillator  potential  HMMUfli 
a  variational  calculation  in  the  thcracftjyneaic  free  energy.    Ills  results 
I  the  thermodynamic  properties  of  a  severely  enharmonic  lett  .th 
lostje  aero  point  energy*    Amplication  of  hie  results  in  v.- 
hindered  by  the  necessity  of  Mflrlnc  Del>yc-type  appropriations.     QMMt 
syaoetry  l         .;  cannot  be  obtained             -o  nodal. 

(h)    Su&aer3  has  also  considered  the  detailed  Shape  of  the 


2D.  J.  JJooton,  Phil.  mc.  }£,  te2  (1955) »  Phil.  SSafi.  ttS,  433  (1955); 


Phil.  Stag.  j«£j  &&  (1,955)* 

3I.  v.  or,  Proc.  Hrys.  Goc.  (London)  jQ,  rM  i 


• 
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potential  veil  ttxaa  the  lattice  sira.    Us  espan&c  tJje  potential  ribo 
lAttice  point  (0,  0,  0)  as  follows : 

w  *  \?0  *  w?(;:2  *  ye  •»•  a2)  ♦  ^jfjF  *  y*  +  a*)  ♦  fe(  ^  i  a  ), 

thus  i$aoria{j  tonas  like    mw»    in  order  to  l^oop  tbe  pro-ale©  separable. 
la  this  egression. 


«  2/72  B  [^(r)  ♦  k/r  f*H*)] 


lb  «  1/2K3C         rI(r)  ♦  6/r  >v(r)] 
°  r 

uhero    X?)    io  the  two~partielo  potential.    S3b©  problca  can  bo  solved  as 
a  etg^crpooitian  of  three  oneniiiaansional  anhanxxilc  oscillators »  rasing 
W0,  ^    and    UV 

SueUer  obtains  an  equation  of  state 


H 


v&ero    E^    represents  the  aero-point  energy,  and  is  a  cocggnUcated  function 
of    V0    and    \l^.    to  ZucUer  points  out,  his  result  is  uneotiafactoryj  it 
is  "baaed  on  an  Einstein  nodel  end  t&erv  G&oroo  the  effects  of  nci£&- 

bors  other  than  through  the  idealized  potential.    Mi  consideration  is 
in  addition  to  the  separability  eggpox&aBtlGn  aade  above. 

(c)    Ducdale  and  rJacXtoaald    calculate  the  energy    E(r)    of  a  pair 


hJ.  B.  Duodnlc  and  0,  K«  C.  ItocStoaald,     Ffcil.  Ua£.  *£,  oil  (IS'iA). 


■ 


of  adjacent  atom  sepere  I  I  rare  gao  solid.    Eb«y 

eoGtne  a  lannarcUJaoae  type  yell  I  MMUVfj      Ltll  ©  fifftft-o 

for  osfoarnonieity  effects,    Saey  MNMtt  that  the  ;;-y 

-bution  to  the  total  enei'gy  ;jortion&l  to  the  Debyc  tera&x  o. 

!2se  Debye  tetapera-;  .;  assumed,  !•  jbe  prop'-  il  to 

&r)  '   .    Using  these  assumptions,  they  obtain  good  ptfsdl;-. 

pNMWM  vs  volisae  i^elationships  sad  cohesive  e:  h  IS*  heavy  rare 

4 

,  and  order-of-naeaitudt  o  * 

ttiey  find  that 

■  *  v  *'3 

g.tt,,». ,    a  F 


M  ••  nolo?  necs,  srons/ciole 
VQ  «  nolo?  voluae,  atr/mrim 
^  •»  molar  enesijy,  cal  /■ . 


2a    r 


jj  -— g  *  K  (Dabye  tesperature,  i  «  lattice 

*  ■■■*         / 


parerseter  ■  ) . 

She  value    L.Ik  :■:  wT^    Batches  eaperiaental  values  for    :<e, 
A,  He    and    15©  ,    with  a  oaxiKura  error  ef  8  - 

Uir  >toa,  and  D^Qdele  and  Msofisn&ld    sb 

clearly  that  theisml  l         "tione  cannot  be  treated  "by  stv&ight 

■  Qebye  tl  is  heats ,  and  that  the  latti.ee  opt:  bflttt 

at  lov  temperatures  BttSt  COM  fron  an.   ■ 

the  dispe..  1  waves  :in  the  solid.     2hei  Je 

th<  of  the  aero-point  notion,  but  they  do  not  Include  a  de.~. 


o:u  cere  then  one  aton  in  tfce  lattice  tyamctes  (other  than  through  collec- 
tive mtion  In  Bound  waves).     In  other  trords,  -U«jy  do  not  include  the 

effect  or  notions  of  the  nei^ibore  of  a  particular  atoa  on  the  potential 
that  aton  that  arioes  frota  the  nei(£ibors. 

Sfoeorlcs  of  Solid  Itellup  Hhree 
She  theory  of  solid  ar  has  seen  considered  Ty  Dornardcs  and 
laa&offf^  (hereafter  tcroed  B.P.)  and  taor  aemardes. 

(a)    3ornardee  and  Prioefcoff  (Doitler-Iondon  approach) 

In  approachlne  the  solution  to  (l),  B.?.  used  the  potential 
(2),  and  a  variational  -wave  function  consisting  of  a  product  of  Gaussian 

functioao  qajB&md  at  lattice  points; 

a 

jL  «  «  const  e  mm  (L(v!)«  (3) 

1  4     '    *> 

If  ^  is  used  to  oval.;    %W  energy  expectation  value,  the 
singularity  in  the  potential  function  yiolda  an  infinite  result,  lb 
avoid  this  difficulty,  B.P.  aodify  the  potential  ty  a  factor  ->&ieh  Gives 

the  Unit  aero  as  r^ ,  goes  to  aero.  Ihedr  effective  potential  is 

*a  ID 

c  a 


WV  "  *• 


e    ^ 


She  cut-off  factor    e::p(-  2  a^/r.T*}    in  chosen  for  ccqpu\  «1 

convenience.    She  value  of         is  estimated  ty  calculating  l;noua  properties 
of  Us    us  /,    vith  various  values  of    0,    until  the  best  fit 


JIU  Beraardce  and  11.  Pstoskoff,  Phys.  Jcv.  13p#  ;}68  (i960). 

°n.  ......    *3v«  130,  1927  {i960). 
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obtained.     In  particular,  the  cohesive  energy  and  i\*3t-aeen-square 

k 

derivation  of  atans  froo  lattice  cites  for  lie  were  calculated  by  a  varia- 
tional calculation  using  the  poranctor  a.  The  valine  of  Q  that  ©eve  the 
noat  consistent  results        -        ■    .     -■.  van  then  assumed  to  apply 
equally  well  to  lie  .  Using  this  assucrption,  they  obtain  the  expectation 


•a 

value  of  ( 1)  for  the  Ifer  saase : 

K(V,  a)  «  H  /  d3  r1  ^*Cr') 


<t(P) 


♦  I  ^  //  «3  ?i  *3  3  **i>  nip 

■  He  WQ  ♦  He  UQ  (10) 

idicre    r  *    is  the  coordinate  of  particle  1  in  n  franc  with    & ,    the 
position  of  the  lattice  point    t§    as  origin,  and  the  into,  a  is  per- 

foraed  over  all  space,    do  specific  voliiae    V    enters  throng* 
She  assumption  has  been  tnode  that  the  effect  of  IfttHtttoalag  syaaetry  into 
the  wave  function  has  a  snail  effect  on    I    and  its  variations.    B.P.  use 
the  condition 

JL  ?{  *  0 


to  find    .  . 

Having  found  a  value  for    a    as  a  function  of  volisae,  B.P.  4m 
TXietrized  wave  function  and  calculate  8  ;anee  • 


J  «  //  d3  $%  d3  »g  f(i2  -  fijg)  ^(r^) 


Uerc    li,*    is  **»  vector  between  neare;;  hbor  lattice  points,    Jy    is 

the  kinetic  energy  part  of  the  integral,  and    Ju    is  the  potential  energy 
part. 

B.P.  find  that    JL«  3.2  x  30"3c    end    J.    =  -  M  ■  -.o"3g,    with 
J  =  -  |  ,  10~3c   »  -  G.02p  K. 

Proa  their  values  c£    E    and    J,    3*?-  then  proceed  to  obtain 
values  for  specific  beat,  entropy,  I  lity,  cad  other  ttaoraodyner. 

properties  of  the  sol 

Shere  ere  no  satisfying  supporting  physical  arguaeats  for  the 
effective  potential  used  "by  B.P. ,  and  for  their  assusption  that  the  ob- 
served  behav  ife    notches  that  of  Be    (through  the  eons  o    in 

the  e.  m  potential  factor).    She  nultlplying  factor  in  the  potential 

energy    ol®{o      c  r"     )    is  designed  to  circuavent  this  difficulty.    Ins* 
ever,  correlation  effects  outside    &«£**    although  ibajporteat,  arc  s- 

eluded.    Other  questions  arise  i&lch.  are  open  to  salon: 

(1)    She  product  of  single-particle  oreitalo  used  by  B.P.  does  not 

include  portiole  cor.  ms  properly.    2ila  fault  of  the  single-particle 

? 

function  has  been  dir,  :  by  many  authors.      Ihe  dirfic-.aty  can  be  un 

stood  by  considering 


i  tfesbet,  -.eve.  Modem  Ffcys.  %%,  20  (Im5i)  dir:  a  the 

relation  problor.  a  aany  references. 


<  flf>  <■  /  $*  (1)  j(p«)    ^(r^)  ^(1)  ^(fite,  4f2  .  (2S) 


If  the  single-particle  orbitala  have  any  ftgUMi  of  overlap  at  all,  arid 
Vi2^r32^    Jmo  a  8incu&n**ty  stronger  then  X/r  ,  then  <  Vm  >  is  infinite, 
contrary  to  £fcysical  actuality*    Shis  consideration  eheds  dotibt  pnrtieularily 
on  the  confuted  exshaaoe  integral*  which  is  detoxrainea  by  the  behavior  of 
the  wave  t  n  in  regions  of  overlap  ♦ 

(2)  3.P.  porforned  their  nstirrtfNrMimi  MMtafl  e.  face-centered 
cubic  lathee.    2hey  nelas  the  conjecture  that  their  maaerical  values  ap- 
ply at  leant  eeniquaatitatively  to  the  bc^~centered  cof  oa. 
However,  since  the  eMchflngB  ^integral  io  hi^iy  sensitive  to  particle  sep- 
aration, their  calculation  of  exchange  la  nuch  leas  li            :>  be  ey^-  - 
cable  to  "both  lattice  structures  at  once  then  their  ca3^ulation 
cohesive  c 

(3)  She  fom  of  spin-lattice  «ttf  eatiferronDcaetie 

crystal  is  sensitive  to  the  e;.cIian{3P  Integral  with  nearest  and  ae:;t  near- 

0 

est  neichbors.      In  a  bc4y-eentereu  t  nearest 

neiehbom  are  only  Igf  farther  removed  than  the  nearest  neichbors,  and 
the  e^chencpe  energy  with  both  nearest  and  no;:t  nearest  neighbors  can  be 
of  the  ooae  order  of  nacnitude.    At  tecs?eratures  high  compared  to  the  floel 
teogporature,  the  ousceptibiiit;>-  t  the  Onglj  iWrtH  lew, 

»  C/(rf  *  Tc).    She  Curie  tcqpe:.  m  by  the  energy  change 

in  a  spin  flip,  including  internctione  with  all  nei#ibors»     (In  Cbe&tc. 
there  is  further  discussion  of  accaetie  offects,) 


UJ«  II.  VenViecfc.     J.  phyo.  radltsa  1ft  U>51). 


)    Bernardes  (Boa-overlapping  *fcse  Function  approach) 

.10  ^  k 

Dcranrfloe  '      has  considered  He*'  end  He    fron  the  point  of 

-v  of  aon-ove^Ia-pping  single-particle  or  ...    9bfl  hosen 

are  square-well  oibitals,  and  Bernardes  calculates  the  integral 

V       -         '  /  *J(i)  0j(j)    V^Cl^  -  ififiJLt)  4f3)tox  *j  *    (S3) 


ucioc 


i        J/l    oin(jr  a  5i  it) 


(14) 


He  obtains  a  sc  ■    (a/ac)    f       ^ '  ^>         <;f)>     *     \v)>   ,    and  win 

aisiag    ^1  ^>  ;     stains  a  value  for    a.     Ike  series  eagansioa  coc  eats 

are  ftUtatttt   tar  "body-centered  MMt  and  »* 

Hence  Bernardes  can  declare  r  ffiost  stable.     He  finds 

Ik  3 

that  both  lie    and  He    asaixae  closest  poe&iag,  wljereao  ejsperUaeatal3.y  at"* 

k 
is  body-cen^cved  cubic,  and  He    Is  bs&a@cnal  close-packed- 

$be  integers!  (13),  aitfaou$3  it  is  finite,  still  overestisBtes  the 

potential  energy.     'Zba  tne  wave  friction  gives  a  different  yueUrtjIlity 

density  describing  particle  ::>r  each  ossuraed  position  of  pc 

•She  sia-ple  p  .  wave  function  Gives  a  distribution  function  for  pc 


3,  Payc.  .12,  J$fr  (195$). 

n.  Be  60). 


■ 


V#  in  c  fossa  »  aRpaofence  on  ttie  pooitlf..- 

SSxLo  dhortconiac  pKiWWIi  an  accurate  eoticsate  of  cvoxveo  potential 

enes&r,  except  in  the  Halt  of  large  separation. 


CHAF2EU  III 

gflS  PIXTOIMS3  KffiHCP 
Meoalni:  of  Particle  Correlations 

u*m   mi  u.n     i  iJt.<i.i>wyiin)n>i>iiiiiw»mMii  ■»■  hi  ■  i»Wi  mwnniiiipni  "   ■*■■■■■■—« 

Aa  eleaentory  discussion  will  clarify  the  neoning  of  the  tozm 
"particle  correlations/' 

She  usual  o.         ;tion  of  a  aany-body  ays  tec  in  coordinate  space 
a  product  0  je-particle  orbitols  (heree  >P«0*)# 

J*    I    ^(r,).  (1) 

Here,  the  fom  of  the  $*o  is  determined  by  the  boundary  conditions  imposed 
on  the  systera  and  by  the  SchrtkUaoer  eo^tlon. 

la  (l),  mo  have  ignored  syametry  requireeientn.     Hi  Merely  state 
that  correlations  reeull  -        rota  the  Paull  ei&lusion  principle  are  soft 
our  concern  in  this  section.      a»  product  (1)  io  the  va\-.  i  of  the 

Zf&rtree  ajpror-lantion  and  io  the  wave  function  of  B.P.~  (apart  fron  syct- 
aetrizations).     *his  type  of  function  *?  erectly  t\ 

state  of    !1    noa-intcractinc  p-         lee*  end  oea  'be  used  no  a  va 


uto  be  precise,  this  prod.  eft  to  Give  a  -....-■ 

eywaetry  under  interchange  of  particle  coo.  fiie  vcnulv,  v 

q  are  unofeeegtdi  -»  *•»  Include  .  , ;    J  ~:  C  P  $,  (a,)> 

■where    P    pezssut**  the      i     ;ea  of  ■etalgnn  the  prope:  or  -) 

to  the  pc  •  t,  qrt  noiiaaliziea     *  the  product.     She  mwrlttMHtj 

any,   U  -uded  I  'Shut  syiabo.  ..JJ.be  consistently  ~o 

■    coord 

a 

--.-dee  and  II.  Priinahoff,  rtr/s.  flev.  11  3). 


.. 
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©tarting  point  for  interacting  particles.     2a  general,  (1)  flail©  for 
strongly  interacting,  deaae  systems. 

Lee,  Huang,  and  Yang,3  in  their  analysis  of  ftam  MBMfttaE 

(hnr&"»sphore)  particles,  use  (l)  in  a  particular  fona: 


»       £ 


Here  the  U»(r)  ore  a  eonplete  orthonoraaX  set  of  wave  functions.  She 
coefficients  0  carry  the  syooetry  of  the  imve  function.  Inis  function 
can  toe  need  to  calculate  the  hard-core  problem,  if  xway  different  coobina- 
tiono  of  It.  are  used.  Lee>  Huang,  and  Yang  transfer  (ft)  to  a  represen- 
tation expressed  by  occupation  numbers      apace.  3ney  perform  a 
variational  calculation,  valid  at  low  densities.  In  a  crystal  order  the 
solution  in  k  space  xaust  transfom  to  a  solution  in  r  space  shoving 
crystalline  order.  Ihie  operation  has  not  been  successfully  performed. 
Besides  failing  to  dooonatrate  crystalline  order,  the  theory  of  lee, 
Huang,  and  tans  is  valid  for  dilute  systeras  only  end  Invalid  for  system 
as  dense  as  the  one  under  consideration  here. 

2b  help  describe  a  state  showinc  crystal  sycwetry  \x>  keep  the 
basic  description  in  coordinate  space  Wt  use  a  less  drastic  separation 
of  $  than  (7)  and  include  a  functional  dependence  on  the  nutual  separa- 
tion of  the  particles, 

!•  ■  ti/VV'  (3) 


3   i  3.  Leo,  k%  Xiuaag,  and  C.  H.  Sang,  Phys.  :tev.   106,  U3  7). 
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V*  study  equation  (3)  W  considering  to  soae  detail  the  behavior  of  tvo 
interacting  particles  in  a  box. 

Bje  usual  Interpretation  of  the  probability  density  corresponds 
to  a  contour  aap,  vith  the  fre^guency  of  observation  of  a  particle  ao  the 
elevation*  and  the  position  of  the  observation  es  the  location  on  the  asp. 
Ttm  "asp"  should,  of  course,  shov  conditions  in  a  3-diraenalonal  solid, 
rather  then  conditions  on  a  plane. 

In  Pic.  1(a),  \m  shov,  in  a  quel  m  presentation,  the  probe* 

bility  density  for  two  particles  contained  in  a  sjtjerical  veil  of  infinite 
potential  and  radius    ra-    She  curves  shov  the  "elevations"  alone  s  iioc 
drawn  through  the  "asp-**    2he  state  is  represented  by    ^(r^)  <{i.,{rc). 
!Ihe  configuration  is  intended  to  represent  the  state  of  loveat  energy 
(©round  state)  in  the  S.P.O.  ep^rasdLaation  obtained,  for  cxipXo,  by  the 
Hftrtree-Focit  aethod.     If  Fig.  1(a)  has  validity,  it  aunt  correspond  to 
vcakly  interacting  particles.    I5e  have  already  seen,  in  the  case  of  3.P., 
that,  if  ve  calculate  a  potontial  energy  value  using  a  9. P.O.  vave  func- 
tion, difficulty  can  arise.     In  the  integral 

<  v  >  -  /  /  l^)!2!*^  •  ?r,)l:V(  |*y  )e?^\  , 

V« 

ve  assuae  V  ■  AMim   ■    -v    ti>  2,  the  integral  diverges  if  ^  and 
$0  overlap  eny**«ro.  In  the  present  case  there  is  complete  overlap,  and 

n  *  12,  so  that  Fig  1(a)  cannot  represent  heliiaa  aton 

We  next  iniprove  the  description  by  retaining  a  general  fujictional 
dependence  on  the  positions  of  both  particles,  vriting  f  *  ^{r  5» 
#(X,  2).  We  cannot  separate  iAl,  2)    into  a  product  of  ^(r^J^Cr.). 


ana 

■ 
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Two  particle  density 
S.P.O.   approximation 


Correlated  single  particle 
density,    f(l') 


t(r{,   r2)lS 


Correlated  two  particle 
density,  particle  1  at  r' 


Correlated  two  particle 

density,  particle  1  at 


Correlated  two  particle 
density,  particle 

Figure  1.   Particle  Correlations 


density,  particle  1  at  r ' 


• 


In  Figs.  1(b),  1(c),  1(d),  and  1(e),  using  this  function,  we  show  bow 
the  'contour  nag?1'  concept  nuat  be  Hodified.    To  draw  the  curveo,  vc  iiust 
think  of  a  whole  aeries  of  saps,  whore  each  asp  corresponds  to  a  part 
ulcer  position  of  one  of  the  particles,  and  dioplayo  a  distribution 
function  for  the  other  particle.     (Wo  think  of  the  figureo  representing 
a  oyoten  with  toe  exponent    a  >  2,    weekly  attractive  outside    t±.  ■  rQ.) 

2a  Fig.  1(a),  ve  ohov  the  probability  of  finding  particle  ((l)) 
[hereafter  labeled  oii^dy  ((1))],  at    $&§  [hereafter  labeled  (!')),  re- 
gardless of  the  position  of  ((3)).    *t)e  eoploy  the  Dirac    5    function  to 
locate  (U))  in    d3rx    at    r\[  ,    una  write    /|*(1,  Dj*  Kl£  •  r^aV^d3*^ 
•  /Ml',  2)|2d3rt  ■  f(l')/  tie  plot    f(i'). 

In  PIgb.  1(c),  1(d),  and  1(e),  we  illustrate  the  probability  dis- 
tribution of  ((2)),  when  ((!))  ia  legalised  at  (!').    Ibis  ft  jtioa 
is  given  by    /|v(l#  £)|L>  o(rx  -  r^)d3(rp  •  \*{V ,  2}f.    She  pitlnat  at 
(l1)  is  shown  an  a  heavy  vertical  bar  of  heitfit    f(l'),    to  illustrate 
the  relative  probability  of  each  of  the  oops  1(c),  1(d),  and  1(e).    2he 
disaster  chosen  as  aboclsaa  la  the  dienetcr  through  (!').    Sartie.::.o        )) 
is  distributed  with  cylindrical  eyoaetry  about  this  diiBaeter.    Uance,  the 
aasaaliaations  of  Pica.  1(c),  1(d),  and  1(e)  are  related  in 
way* 

In  constructing  the  figures  we  have  aaouaed  a  strongly  repi~ 
interaction  of  radius    rQ.    If  the  Interactions  were  weei;  everywhere,  til 
representation  1(a)  would  introduce  a  snaller  correlation  erro 

For  three  and  store  partlclos,  we  proceed  in  the  sane  fashion,  I 
with  a  new  set  of  naps  for  each  additional  particle.    For  oxoaple. 


three  particles,  Flga.  i(c),  1(d),  end  1(e)  can  represent  the 

tioa  or  ((!))  end  ((2))  for  a  £iven  position  (3')  of  {(3))*    We  can  vrite 

the  three-pa:  function  as    y{l,  2,  3)>  n  three-pa:  an 

.,  8)  $(2,  3)  $(1,  3),  the  tvo-particle  app.  EttOBJ  or  as  ${1)  4? 

$(3),  the  alogl#~portlels  appra;:iaatioii.    In  the  spirit  c  g  » 

inprovermt  to  the  r..P.O«>  and  since  we  vill  have  no  three-par  ■ 
actions  to  deal  vith,  om-  work  usee  the  tvo-particle  epproxlr.; 

1(1,  2,  3)  ■  ♦(!,  2)  t(2,  3)  *U,  3), 


.,  in  general 


id,  2, 3; . . .  a)  -   1   if(i,f  **y 

1  >  $    *    j 


Use  of  the  Interpartlc^  Coordinate 

tpi  Provide.  Conflation  tfc 

■a 

She  cathode  previously  applied  to  the  ground 

vhioh  kt  a  descript  '  the  ground  fttteti    tat  OKewplt.'  ), 

ress  correlatiODA  through  a  produ  .^lc-partici©  orbitals.     As 

described,  in  Chapter  II,  such  a  descripti.-  |  poten- 

tial with  a  stress  singularity  at  tl*e  origin.     Including  the  iaterpa-: 
&  Jao  state  function  would  reaove  this  ol.:.  fet  two- 

particle  function  of  Chapter  n  offers  a  way  of  dolus  this. 

*«t3  considered  this  pa  .In  parallel  caXcul 

tione  of  the  temperature  dependence  c 


4  5 

Price    introduced  a  fixed  interperticle  dleteace,  **u         t    esotsaed  a 

hard-core  repulsion. 

G  7  8 

Jaotrov,     «vlloo,     Jwaaoto,    and  Xtetachaen,'  have  developed  a 

variational  method  (the  Jantrow  Esethod)  noinc  a  trial  function  IT  ^(2^  ) 

TA 

n    f  (r,  -  r,).        Sale  function  reaeiables  the      I!    ^{r,  r)    of  the 

loot  eection.    She  Jaotrov  saethod  includes  a  aerie*  expansion  in  clue  - 

integrala  ela&tar  to  tisoee  used  in  obtainine  tbe  partition  function  in 

11 
classical  statistical  nochanics  of  donee  o^tea*.        Sic  cluster  nethod 

of  obtaining  the  variational  paraoeter  <  la  U         *i  to  tbe  eolid  etr. 

because  it  nesiaee  liquid  disorder,  and  tlserofore  tbe  Jeetrov  aethod  will 

not  serve  our  needs.     However,  tbe  method  has  points  c 

our  aetbod,  and  we  will  return  to  it  later. 

Bruedoier     has  developed  a  nethod  ;       ig  two-particle 

iaterectiono  correct fy  to  arbitrary  order.    2be  theory  is  a  eotterlng 


**P.  J*  Price,  Fays,  .lev.  S&  2fy)  (1955). 

50.  K.  Hot!  Fhye.  Bev.  >). 

U  Jaotrov,  Pfcya.  Ucv.  gg,  IfcP9  (2£$5)# 

v.  3.  Aviioo,  Jr.,  Ana.  Pays.  %  251  ( 

°F.  Juodc  .  -  Theoret.  yoto)  1£,  J). 

F.  Iweraoto  and  H,  Xaaada,  Prog.  1iliWi1v<  Hiyo.  (ilyoto)  1£,  5^3  (1957)- 

«  J.  DwwSaaoa,  Pays.  2»v.  121,  6%8  (ly6l). 

2ne  author  io  Indebted  to  Professor  Volte?  Kobe  for  pointing  out 
tbe  connection  between  tbe  present  eeafe  and  tbe  approach  of  6,  7,  0,  ,  ■ 

11 
J.  F«  Meyer  and  II.  Q.  Mayor,  atlcal  Mechanic*  ( ' 

Sons,  Inc.,  Sew  Yc  1)0)  p.  1/77- 


*%•  A.  Botbe,  Up*,  ftpr*  &£,  1353  (38*) 


theory.,  and         -   jq»  an  expaae  Ml  wty 

applyinc  the  iNidtoMtf  oethod  to  a  solid  lice  not  been  developed.     She 

2 
pseudopotential  asthod    chains  this  aieadvantaee. 

is 
PI.         -o,  J  In  studying  the  atactic  apectruv        h<l;lM»t  used  the 

set  of  coordinr.  cd  ear  .lie  Pluvinaee  added  the 

feature  of  including  the  iaterperticle  AU    .        WjUn 

order  HMd&tOttftAfe     Bu  acro-or&cr  vavo  function  is  thus  a 

the  space  coordinates  of  the  particles,  conta 

tance,  and  cannot  Ixs  separated  into  a  product  Of  le-partlc.'  >aa. 

It  possessec  none  of  the  disadvantages  ascribed  to  other  aethoda  in  the 

:cdiug  para^^aphs. 

Uoiiit;  Piuvinace'c  oethod,  Walsh  and  Borwits      (berec 


y.B.)  have  tronalfli  be  anwjp  levels  of  e<.  atons  and,  Sn  t 

yet  better  results  in  aero  order  ti»an    .: 

calciilatioae  usually  piiftaajwd  with  hydros  ,  is.     The  method 

B  also  be  applied  to  our  problm,  vi\.  -.Tonaale  e:pectation  wt  equally 

good  results* 

She  Pluvlnagse  Method 
2mj  Piusviaocse  raethod  introduces  the  interpartlclo  distance  into 
the  nany-body  lie.  gas,  mi  •  tt0tlOA#  yielding  an  exact  foroal 

*^«  Piuvinaee,  tab  Hmu  2*  I2*?  (1v50). 

ts,  Pfcys.  itev.  U£.  1206  (1959);  Bsys. 

z*r.  yg,  127^  (i960). 
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esproooion  uaicn  iacludee  *3®3Aolt  (Sependence  on  all  the  intxaportelcXo 
dlstancee*  fto  -noil  an  She  opatial  pcssitiono  of  all.  the  pc         ■*« 
Uo  vslte  tbe  amv-bcx^  lin£ltoiilan: 

3e«j  ve  describe    n    Identical  particles  of  aoso    a,    interacting  wi" 
tbe  potential    V(r,.  A).    Vie  have  alreaty  seen  that  ve  can  take    V(r.  „)    to 
be  tbt  lenaerd-Jbneo  potential 

V(r^)  *  ^((cr/r)12  -  (e/r^J 

f  m  A  (5) 

V(x% )    is  an  J         ially  applied  ffllmtlalj  the  se&»  ftraction  of  space 
for  all  part-aleo. 

tt»  vill  nov  follow  Pluvlnci{3e  la  introdusias  the  interpartlele  dis- 
tance Into  the  Ifcrailtonian. 

In  any  calculation,  law  a    independent  variables    ( u, ,  v,, 

**j  y^/  vo/  w0,  .».)>  «e  can  introduce  another  set  of    &    variables    (:*_, 
y1#  a.,  »0,  I  ...),  related  to  the    n     by    x    »  3^(u  )t  and  pooaibly 

simplify  the  calculation,    tfo  can  even  have     m  >  n,    provided  that  vo 
aloo  provide    n  -  n    independent  constroi:  the  function  u  ). 

In  f®pmt&&  A  m  ahov  that  the  following  statenseat  Is  true: 

If    ^jCiT)  »  ^(n-,  v,,  v,,  \Vj>   ...  vj),  i  function  c  .  tesian 

>idlnates    {**$  v, ,  v.)    of    n    particles  in  ft  cyotcsa,  \m  &ny  set  <,$> 


..VI' 
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onothcr  function    fgC*  )  *  tg(*j»  *j.#  si  ***  ytl  aU'  r12  r13  ""  r3-l,' 

which  contains    (U)(II  -  1)/         ■.=»  coordinate©,  and  we  iaay  require  that 


$-(u  )  «  ^(s  )  everywhere  la  space.     In  orSor  to  do  this,  we  must  have 
(3*1  •  (H)(U  -  l)/2)  iratogpoacient  relations    x1  *  x^u1),    ana"  the  rule 

jtfffrHNl41 

i?    53    i-i  i? 

In  our  cose,  the  8o?^ltoaiea  contains  the  iaterperticle  distances 
r„ .,    so  we  -ore  led  to  latarotae  thex*  as  c  x&aatea. 

ifc  then  hove 

Hm\    yi"vi     ai°vi 

(6) 

rl£>  "^(uL  -  ^  *  (vl  "  V^  *  (V1  -  V 

as  the    38  •  H(IJ  -  l)/2    ladgpeadeat  relations* 

la  what  follows,  we  Shall  us©  c  vove  function  of  the  type 

which  gives  the  seao  results  as    $, ,    simply  because  we  roq 

tfo  chall  find  that    yn    yields  a  useful  natural  separation. 

%  drop  the  distinction  between    x-     sad    u~,    hut  ramitjer  the 

special  role  the    k*t    play  In  the  functional  flora  of    iu*«    Fran  fggorn 

A  we  heve  the  rule  of  differs;  a; 

4  -  4 +  i  &  ^  <7) 

liero,  to  the  ri$it  of  the  arrow,  we  speak  of  a    ^0    type  of  doc  m, 

to  the  left,    $.« 


..    I 


Front  Append:  N  neve 


■;,     -» 


. 


-  - »  -  ip*         | 

L 


^    3    •     -    I 

■  ■>«■  ■■«■  ^Mjfcw^^m  i»fcnfrw«nffi      »w*iw» 


*  s 


-fj. 


*  °i<3 


U 


U 


Berc,  we  have  put 


4J     3o„  „ , 


s,ij 


I 


MJ 


«     (::,  •"j^), 


■fij 


(10) 


Here    e*.    is  the  vector    r<  -  .:•. .     it  is  uaed  to  tmpwii  relative 
placee»nts  in  tenia  of  the  variableo    ;:..   ■  u.  >  %.  **  *«j  ft  -     S&e 

eyrabol    V4 ,     Ui  one  0*  the    n(U  *  i)/2    a&ait.  ;oardincteo, 


bl  •  il)     •  (v.  -  vi)"~  ♦  (v    -  ILr  *     Hi  ^e©a  b-  -    end    «t, 

since  JMfil  Itll  1  W*l  in    x    space,  lib©    bfim    or  ,  ,  operas  ;e* 

.a 
tiona  of    s,  y,  ft*  cad  o    but  not  on  fraction*  of    r,  ..     likewise,  derlv- 

.ves  ouch  as    d/3r4 ,    do  not  operate  on  :.ono  of    j  do  operate 

on  fuacti-  ,.     In  addition  to  the  o  vector    a,  4*    ve 

ij  iO 

-will  aleo  have  need  of  tbe  pc  ■  MP  no.se  of  pB  j, 

ttqpoMHUMd  |  f$  1    space* 

(11) 


->  «*  4* 


We  can  write  the  Iteslltonian    II, 

II    *   Uq     ♦ 
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"o 


.i 


2 


r  i 


H' 


.  t£  !mJJe*m  a 


a 


kiJ 


• 


U 


2  s2 

^  £    8lj  '  *ifc  V^r^ 


ilfc) 


1X1(13),    ^- 


n< 


~j 


fU   ^ru 


An  atteapt  to  use  the  total  asailtonian  to  ctetenaixie  a  set  of 
ei<$enstates  of  the  Hchrodin^ar  equation  appears  uninviting;.     Instead,  we 
use  em  approxirmtien  introduced  oy  Pluvinage.    tie  look  up;  as  a 

perturbation,  and  look  for  the  o totes  of    ft.: 


Sere  |  is  a  function  of  the  3B  ♦  8(!f  -  l)/2 


*!'  *2*  *** 


r»*  rl3'  •••  r23*  rtM  ***  r3U'  iM  ri»l  I  '      ^  •*tiUk  <*5)  by  the 

standard  raethod  of  separation  of  variables;  the  form  of    VL    Bakes  the 

solution  separable  as  follovn 


S  P 
P 


J  L  *  ^  j       — 


(16) 


Here  ve  show  the  spin  function  *j.  f  a  function  of  the  spin  variable 
C.   Ihe  state  function  <§  It  an  approximation  to  the  true-state 

tion*  If  the  perturbation  potential  H'  is  truly  ooall,  the  epproxusa» 
tion  vlll  give  little  error  In  the  eaergv  ievei  of  the  ground  state. 


i*    'J«U 
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In  vor  ith  the  function    J  ,    ve  auot  remember  tlje  role  of 

poraoeters  of  integration,    We  oust,  therefore,  follow  a  program  vhich 
lias  the  following  oteps: 

(a)    tfolve  (15)  for    J 

(a)  la  using  ^  and  !!',  let  ^x  operate  on  x,  n  ,  t 
only,  end  let  ^/^r4 ,  operate  oa  r  ,  only. 

(c)  After  differentiation  [step  (©)!*  introduce  the  constraint 

(d)  After  Introducing  the  constraint,  evaluate  nx  timet* 
"ay  into£ratian. 

Uote  that  the  order  of  steps  (b),  (c),  end  (d)  cannot  be  altered.  As  in 
any  analysis  involving  a  pareneter,  we  con,  however,  firat  introduce  the 
constraint,  and  then  perfom  any  itraalalog  differentiation  by  ordinary 
rules  of  differentiation  not  Involving  &/ta  ,.  ite  shall  find  that  this 
latter  progren  will  perait  the  calculation  of  the  average  cohesive 

Pluviaaoe  and  W.B.  have  shown  that  (1?)  yields  excellent  values 
of  £  for  several  systeras  they  have  analysed.  Further,  they  have  suc- 
ceeded in  obtaininc  (pod  resulto  by  conslderine  oaly  EL,  after  oniiing 
additional  assumptions  eaoceroine  the  function 


ao  r> 


i*  i   . 


She  Haxsiltoniari  ill)  still  contains  the  singular  interpnrticle 
potential.    2he  various  applications  of  the  Pluvinace  nethod  differ  in 
the  handling  of  the  potential    V(p4  A    jfcro,  we  ahov  hov  to  use  the 
Pluviaage  aethod  vith  strongly  repulsive  fl 

&    In  (15)  con^n.  a.  ftdtafec  t«»  aw*  on    V 


S  P    c 

p    i,j  L 


1 4u  •  i  W 


tie  can  natee  this  tern  vanish,  and  vith  it,  the  potential    V4  -(*«*)# 
res^lrin^ 


£*j*i«4j> 


X(riJ)s*0'  (lB) 


Expression  (13)  can  be  looted  upon  as  a  nathenaticol  art:  e- 

slened  to  sumouat  the  difficulty  of  the  singular  potential.     It  MB  aloo 
be  considered  a  Scbrodineer  equation  defiolnc  a  particular  acc&e  the 

coutinuua,  or  states,  the  lowest  iyinc  of  such  states,  that  of  aero  energy. 

In  Append!::  B  ve  solve  (13)  by  changing  ftm  variable    r    to 
variable    q  *  c/r,  tatting      cc  c  power  series  in  ascending  positive  MM 


^fe  could  write  the  ri<  .;ide  of  (l£)  as    e    .  *-  i^V*)* 

Wtt  would  find  no  bound  states,  hence  ).    For    e   .  >  0,  there 

would  be  oscillations  in    X    at  lar^e  distances,  vhich  would  cause  real 
trouble  in  des  ;  an  extensive  solid.    S&aoa  wo  hope  to  have    . 

tional  faro  of      (r,  ^)  be  independent  of    r±.,    ve  oust  tlftHtH  i- 

latioas  at  the  outset  by  having    €±.  =  0.    ifo  could  still  solve  (lo)  vith 
*U*U    aa  ^  r^%t>  ^^  *ajen  9et    6  «  0,  but  this  process  would  on3y 

serve  to  -.luce  additional  teras  in  the  series  expansion  for  . .  ii 

would  contain    <s    on  a  auitipiier,  a  needless  cceplication. 


of    q.    She  result,  for  the  physical  pareewters  of    Xle^,    ie  a  rapidly 
converging  series.    The  two  linearly  independent  solutions  4 

0  "  C0  F0(flr/r) 

^  «  C^o/r)  F^c/r) 

F0(^)  •  1  .  Aa*  ♦  sP/(3  •  7  •  *  *  3)  x6  ♦  B/C1Q  •  ,)  a:30 

-  D?/(il  ♦  12  •  0  •  7  •  fc  •  3)  ■?" [  ♦  .  -  . 
rx(»)  -  1  .  S)/(fc  •  5)  ^  •  2?/<9  •  S  ♦  f  ♦  'U  £  ♦  B/(U  •  10)  y?°  ♦ 

Hero    D  *  n  a*e  X  M*  *    which  has  the  value  l£.60  for  He   . 

At  large  r,  *vf*cn>  »*nd.  v4sea  jr-*a  froa  above,  X.  tends 
xwpldly  to  aero.  Shis  behavior  ll  characteristic  or  a  hard  core  interac- 
tion. At  large  r,  "X,  •  0;  it  then  M  a  Marfan  at  an  intexaaOLetft 
value  of  r,  and  hao  a  finite  valuer  at  r  *  o\  In  Appendix  3,  X,  is 
shown  to  yield  a  oolution  vith  an  outgoing  at  large  distances.  Ebls 
is  an  unsatisfactory  chexocterlotic  nahich  cannot  be  penaitted  in  jL.  VJe 
thus  choooe    C,  *  0. 

We  now  dssaoasfcrate  that    CQ  *  1.    Since  the  Eeaiitouian  cloca  not 
contain  the  tiae,  -we  aay  choose    X  to  be  real.    Cinco    f    is  real,    CQ 


oust  be  real. 


C0  "  °0*' 


At  large  fiietances  of  separation  the  interaction  has  no  effect  on 

the  t'wo-particle  nave  function.    Therefore,  the  protafeility  den:- 


•37- 


particies  i  and  J  1b  independent  of  r.  4,  and  raultiptying  the  wave 
function  by  the  equate  of  the  correiation  amplitude  leaves  the  wave  func- 
tion unchanged,  dat  is, 

Q  O 

therefore,  C^  ^(0)  »  1,  or  Cq^  »  1.  •!  taJse  the  positive  root  of 

If  we  were  to  use  the  \&B  approcdiaation  to  obtain  the  wave  twn 
tlon  for  r  <  o,  ve  would  find  that  X  would  continue  to  fall  c 
rapidly  oo 


esp)  -  /    CV(a  -  r)2a]   h  d  or 


x  [V(a  -  r)2a] 


Hence  ve  ore  justified  in  assuming  tfcat  is  vanishlngly  snail  inside 

r  »  a  (provided  oni^r  that  we  use  caution  in  calcinating  derivatives  t 
r    S  s).     our  representation  for    /,     the  aeries    F0(a/r),    has  poor  con- 
vergence properties  for    r  <  a.    Because  of  the  ssaallneas  of    x(r)    for 
r  <  o,    we  can  cut  off    fAq/t)    \y  writing 

Ar13  «  F0(B,  r^)  9^  -  rQ»  (») 

Here  rQ  £  a    where  rQ    cannot  be  ouch  less  than  a,    and  0(x)  is  the 
step  function  0{x)  •  1,  x  >  0;  9{x)  *  0,  I  <  0. 

In  Fig.  2,    *'*£.)>  *(**.)  j  oa*  X^ri^X^rij^rii  ^w*  P1**"***2- 

■a 
for  Her.  In  Pig.  3  ve  illustrate  the  ccqparison  between  the  series  func- 

l  FQ   and  two  appro:     aii&Jytical  representations.  'Hie  function 


mivLnm 


r      -r(r-r0) 
1  -  — -  e  shown  in  both  Figa.  8  and  3  lo  an  extremity  cIodc 

Shis  Id  the  trial  taflft&tt  v&ich  hois  been  fovniQ.  to  give  the  lowest  energy* 

6 
out  of  oeay  tried  in  the  Jcstr.  uxU 
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CEAP3EH    IV 

MMBr-Boay  picaiEH  m  spued  ieijih  mag 

.1      ....III     II. » I.      I       .,1.    !■■,■    ......        I        I  I .11.1, I|.        «■■! 

asoalty  Ustrir,  Poi'nulation  of  the  tlaay-3ody  Wfcws  function 

la  this  chapter  ve  shall  develop  a  tfceory  of  the  rasny-bo&y  problem 

■which  can  be  applied  to  dense  solids,    2he  developnent  uses  the  densi 

1 
natrix  fomalisn  of  Xowdia,    vhich  w*  review  briefly  here. 

%  consider  the  sie&y-body  v*ve  function    t    xdiich  fulfills  the 

syaoetry  condition. 

p  tfx^,  x,  .  .  .  %)  »  ( s)p  *(£  ■  I  t,)  (1) 

vhere    P    io  a  peirautetion  operator  working  oa  the  subscript.: 
ordinates,  sad    p    is  its  parity.    For  aose  particles,    q  *  1.    £or  Fermi 
particles,    r\  ■  •!.     She  spin  variable,  if  any,  io  included  in    X.    A 
physical  quantity      Q    associated  with  the  system  la  represented  by  a 
Hensltlen  operator    0      .     It  aay  be  expressed  as.- 

vhere  we  separate  the  operator  into  its  zero,  one-,  ta&»4  .   .   .  particle 
operators.    She  prloe  oral  to  teroa  having  two  or  aore  equal  indices. 

Tb  evaluate  the  average  value  of    Q      ,    uo  introduce  a  series 
density  natricos  of  various  orders: 


5)  J  Pbyo.  Rev.  2L  W  (J055). 


-In- 


r(i{\\)  -  m>*U'23  •  *  ♦  n)  f(K3  .  .  .  it)  c5x  .  (3) 

Ihe  syntool    tiL     cicaifioo  int«  on  ovor  all  variables  ftsoupt    !•« 

KX*  J^^)  ■  $2^  /  '3  .    ■    . 0  ^(123  ...  3) 


r     <x£  g£  •  .  .  z'\::x  ^  .  .  .  x  ) 


-ir^-rrr/ t#(l*2'  ...  p'  •««  ■>  f(M  ...  P  ...  0(4GL  ...  d£) 

p.vy  *  P^»  -  P 

r^hx^x^   .  .  .  ^LilX-XL  .  ♦  .  tj  «  ^(Wj'  ...  U')  4f(X23  ...  If). 
She  aatriceo  (4)  fulfill  the  following  condition.. 


_i.   .-* 


But      XW>  i[\^  \)  m  r(K|  ^l^ll,),  regaMloss  of    t}.  (5) 

l/e  also  have 

r         (Xj!  ~C  .   .   .  X^l^  Xg  .   .   .  ^J 

~  inHhrp  /  r(p)<%  t—  %*  fPir  S  •  •  •  Vi  %)aV    (6) 

She  diagonal  eleneata  are  or  opc-iol  ixsporta- 

r(^,^)ar(£^y.  (7) 


Shcoe  natricos  are  syiaaetric  for    %  m  ±  1,    Ml  con  interpret    rC^dv^ 
eo    II    tiae*  the  probability  Of  flatting  a  particle  la    dv^    at    *.|    with 
spin    £.,    regardless  of  the  positions  and  spins  of  all  the  other  gM 
tlcles,    Cin1,1arly,    r(&,,  ^dv^  dv2    ls  ^  *»**»?  of  pairo,    5?(n  ~  l/2, 

os  the  probability  of  finding  one  pnrtlela  la  dv,  at  "r,  vith  spin 
£-  end  another  In  dv0  at  1%  -with  spin  £->,  regardless  of  the  spins 
end  positions  of  all  other  particles*    According  to  (3)*  vo  nave 

She  matrix    r  **      is  entis^aoetric  it:  each  act  of  indices  vhea    n  •  -  1. 
therefore,  if  two  or  more  indices  MM  c<r>ml,    I**"  ■  0.     2n  particular 

Ti\,  \)  *  0,    P5'(^f  £>»  l£,)  •  0.  (9) 

xhesc  relationahipo  are  egressions  of  Hm  Null  principle. 

We  now  consider  the  uee  or  operators  to  obtain  average  values  of 
physical  quanti      Usinc;  a  tuo~particie  operator  as  an  e;«ople, 

</  ay  -  /  *#(s'  a^ft  d:tx  .  .  .  4. 

Svaluatine  the  inteeral  for  the  general  pel*    21;,    vo  have 
<£ ft  7»|  n(n  -  1)  /  $*  p..         ^  ,  .  .  1 


/  ftru**3  ...  if)  %,  *(123  •••  n>}  i»««?|  i 


.'■■■' 

.     •     . 


Applying  (3)>  ^o  est 

&     «•    K     V 

y©  use  tfce  con. I  I  tfeftft    QL  ,    operates  only  on  t&o  iuprir.xs&  cooi-dinates 

X.,     end    *.,    end  not  on    X,.%  X:.      Bence  the  noa-fcLaeonel  elaaents  are 
needed  to  er^rttOG  tisc  ercpectetion  values  of  differential  operators.    Doing 
tlio  convention  of  (K»),  we  #et,  for  the  General  operator  (2), 

♦  /  nj^xjf  ■'"■  |x  i  L)^  ax^  dr^  ♦  .  .  .  ,  (u) 

In  our  problem,  ve  neve  the  two-perticio  ope  IK  *  rjl^' 

end  do  not  Introduce  th^DO-pert-v-lc  operators.    Hence,  a  description  of 
the  solid  in  tense  of  tenuity  oe-trlcoo  .Thould  give  occurae* 

values  of  <^*%p>,    oo  follows: 


In  (22),  ve  taiUe  cere  of  uoimiisation  by  t&m  tenoGdnoctor,  unid:  -a- 

•acted  to  0*  unity  to    ^  -  1. 

•tac^-Particle  Density 
(a)    In  this  section,  ve  use  tlao  fremiti  :e  density  oatr 

intixxluced  in  the  praceding  section  to  build  up  en  I?    particle  dose  rip- 
tion  of  tije  col.    •      fflbio  dm.;           Lon  vill  be  8££flr  -that  It  will 


. 


> 


■'^Mli 


~&5- 


dcscribe  the  psropertiea  of  the  solid  in  the  oiaGlc-partiale  sense,  $uat 
as  the  £  later  deteralnent  describee  interoctine  system.     Our  de- 

scription toI%o  the  correlations  explicitly  into  account,  on  addition: 
reflnooent  not  included  In  the  Olatc  tion. 

(b)    Sing^- Particle  RWilfcill%  Dousity.    We  atoall  write  dovn 
the  B*body  density  ijatrii:  in  terns  of  the  full  Pluvinage  -wave  function, 
and  reduce  it  to  a  oinGle-pe..  i  i  .ty  lartdlx.    T,Je  then  show,  for 

dense  systetao,  that  a  sirgOe  representation  con  lie  used  go  the  solution 
to  an  integral  equation,  yielding  a*  for  the  single-particle 

density  oetri;*,  %&ich  includes  correction 

She  B~body  density  antrix 

2a  (13)?  vo  show  the  diGGoaai  forn  and  do  act  introduce  priiae  coordinatee. 
She  result  is  1  ".ently  General  for  our  needs,  since  ve  do  not  seek 

expectation  values  of  differential  operators,  but  txuxjly  the  probability 

density.    IJon-diacjonal  complications  will  concern,  us  in  later  sections. 


When  ve  use  (13)>  ve  trust  reaeober  t  of  const..  r. .  »       ;M*®ii  ■ 


It  io  not  necessary  to  specify  the  forn  of     '    until  ve  attcsgpt  to 
uuaerieal  calculations.    "Me  shall  aasuoo,  however,  that  U  ft  1     versa! 

function  of  the  distance  between  any  two  particles,  independent  of  the 
spin  of  the  particles. 

*toei  Bq.  (3)  t&e  siocle-particle  density  zaatri:;    r(£.. )    ana  be 
obtained  ..ocpntinc  the  n~x»rticle  Mft  or  oil  caordinates  c 


■ 


r(\)  •  (coast)/  r*(S^ 


- 


*  •  •  • 


<*) 


Zio  constant  ic  a  HillllUTHfl  constant  dependinc  only  on  the-  BBMi 
particles*.    £a  the  hapter,  v©  twUmtjj  it  cm    K 

iteeli. 

Wa  can  rewrite  (lU)  in  tanas  of  the  representation  (13)' 

r(i)  -  (conafc)  l  ^U)  ^(D  /♦*(2)  *2<£)  ...  4g(a)  *nUD 


(15) 


2a  (15),  \m  have  explicitly  ghovn  ttic  diaoonal  fora  of  the  tksacity  astr. 
!2ie  poscEutatioaa  refer  to  tnterchatigiag  tbe  cocr&iirateo  or  the    $,. 
tooth    $?    oiid  ^     cinultaneouoly*    Sn  reoultlne  function  is  oymaabr 
the  interchange  of  nay  two  particles.     7bc  arjtiOQl    IL  ;r  ;'or  qynactir: 
ohovc  talc  ajHaetxy.    like  lent  relctionohlp  la  (5)  ohovs  tnio  property. 
An  aoa\3qption  ftl  irspllecl  !  king  (15)  la  t&io  ra$e  la  that  ue 

are  neclc  the  spins  of  the  particles,    tfe  luwe  eaouaed  tbat  If 

4(1)    4(1) 


,(2) 


(l'c'lia)  « 


KUt)  *a 


*<2f) 


V) 


ra 


♦n<a)    *b 


the  i 


x*s,(aa|j2)  -  |*a(i)|2  \fy*)f  *  |*a(a)|2  |*b(x)| 


2|y  co  doing,  vo  hove  dropped  ftm    r    '(32  J 12)    on  overlap  t. 

-  [(♦*(!)  $a<2)  $£(2)  $b(l))  +  (oaragae;;  coaJUGBtse)].     "Enis  tors  contributes 

(2) 

positive  or  negative  ooeuato  to    r      ,    dependlne  on  tfre  relate 


functions  of  $    and  $.  «     »fe  shall  find  that  (a)  correlations  keep  our 
particles  voll  separated,  and  (b)  the  $'o  are  orthogonal  for  particles 
veil  aeparated,  oo  that  contributions  o;      m&*§  ere  very  snail  eciegpaxed 
to  terra©  like     {$..  (i)r .      Hence  ve  ore  .justified  in  ng  the  overlap 

tern  uhen  It  appears  in  (15).     Ml  vlll  i*ind  that  v®  cannot  dr.  -a  tens 

\jfcen  m  consider  e:  whence  effects  * 

All  ascmae  throughout  that  the  WWBtffie  pmmAm  Of 

diagonal  elceaents  ahovn  in  (&S)  io  adequate  to  describe  probability  deao* 
ity  effects,    tic  ohall  find  later  that  exchange  of  recto  are  aocIiG~^# 
thus  Juotifyinc  the  assumption. 

Uriti;.  )  in  aodlfied  fom,  ms  MM! 

Til)  *  (const)  Z  ♦/(!)  /{«  Pr  *22(2)  ...  ^(j)  ...  */(») 

^3  *   *   '  :<B-HI}  ^  *»  ^  ' 

!Ihe  tera  in  curly  brackets  ie    r"       '($,  2,  ...  J  -  1,  3  ♦  1,  ...  II). 
Because  the  diagonal  form  of    P    io  qpMrtSle  on  inter^hasee  of  lt§  co- 
llates, the  ordering  of  the  labels  U  the    $'s    as  veil 
as  the    X'e.    Hence,  in  each  nenbe          olio  Mi  over    J,    the  nultlple 
Integrals  are  identical  in  fom,  independent  of  the  position  of  the  in&«K 
J.    \te  con,.  therefore,  factor  out  the  into.       .    a  each  s-jn  over  J: 

r(D  •  (const)  f  t  <^2<i)  ]  /  r**"1^ j,  e,  . . .  •)  |  K^  v^«    U7> 

tfe  oust  now  consider  th< 
He  assume  that  w  out  a  m  p  %i  the 


«4o» 


nony-partiele  description  (13),  or  that  there  ©rjists  a    r(l)    K  -u- 

scribes  the  average  behavior  of  the  solid  by  aaeicninf,  certain  props.. 
to  a  ain&Le  average  particle  labeled  (l).    !23e  refinataBixt  that  %aa  build 

r(l)    should  include  the  effect  of  goyyleittaqg  between  the  particle 
with  coordinate    r«    and  the  avercQs  position  Of  all  the  other  pa.  j, 

but  will  not  include  correlations  within  other  pairs.    (Shese  other  pairs 
would  enter  into  a  two^body  density  aatri:;  description    Til,  2),  which 
describes  {»)>  correlations  of  1  and  2  vitfa  the  average  positions  of  an 
the  other  particles,  and  (b),  the  actual  correlation  between  1  and  2.) 

The    ^(*\v)    *«*  (*?)  carry  the  ..iafclor.  .-otueen  atoc. 


1  and  the  other  stems    fc*    but  ve  need  to  MM  average  positions  of 

the  other  atom,    For  this  purpose,  expressing    I* '  as  a  jjinjinit 

ainele-partlclo  density  nsntrioes  is  adequate. 
ibercfore,  if  vo  set 

Mi 


FW(1,  2,  3  ...  3)  -  -4      I    rU), 

tr   j»i 


we  httW  en  ecrprenaion  \&ich  satisfies  the  relations  (3),  (h)   and  (5)  ap- 
plying to  density  aotricee,  and  specifies  the  positions  of  the  parti' 
adequately  for  the  purpose  of  evaluating  the     al  in  (!?}. 
Hn  can  now  write  for  (17) 

r(D  *  coast  1  <k2(l)/  1  rit)    V    >:     >  d£ 

(15) 


«  const  |r,  $/*(!) 


U  J     M 


i  /  rdc  • 


. 


Equation  (13)  is  en  intOGrol  equation  in  y(<5)  which  gives  the 
projection  of  a  aingle*p«rticle  state  frcta  the  ouTy-body  state.  'So  ob- 
tain it,  wo  have  ossvaed. 

(1)  aJjo  tvo-partielo  separation  of  the  wave  fgnoticn,  (13):- 
an  adequate  description  of  the    H    body  oyotca; 

(2)  Aaons  the  sinGie-pcxticlc  states,  overlap  haa  a  ifHjIlpH 
«f..'  le-partielo  properties,  in  that  intecrals  of  the  type 
/  ♦j(l)  %(l)^v  k  J*  <3    all  vr. 

(3)  \j(**t)  la  a  uiivorool  fuoct4on  of  the  distance  r»., 
roGarOleeo  of  the  identity  of    i    and    j» 

\te  hove  not  aosvnod  boundary  00Q  M  for  the    $'s,    nor  a 

©pacific  functional  foan  for    X.  4{r.  4). 

Si  the  neirt  section,  uc  shall  use  (l£)  to  Obtain  an  explicit 
single-parti ale  representation  for  solid  lie  . 

A  Golution  of  the  Density  Equation  for  Dense  Bolide 
It  cocao  obvious  that    r(r.)    in  (Ik)  is  a  function  which  has  a 
CTMlrmm  in  its  spatial  dependence  et  1;  points  and  shows  the  trensla- 

tional  eyaraetry  of  the  lattice.     In  fact,  we  can  define  a  solid  by  requlr- 
ine  that  it  be  described  by  a    r(r4)    with  these  properties.  .vc 

additional  asauriptioas  that  to  a  good  order  of  appraxir»atio&    y(r4)    has 
spherical  syanetry  in  its  behavior  in  the  iape&iatc  neighborhood  of  lat- 
tice points  (tiiis  consideration  ie  independent  of  the  fact  that  the 
oaxhaa  of  the    fCrJ    are  disposed  accc  to  lattice  eynaetgy),  that 

r(r .)    decreases  raonotonicaliy  away  Ifem  lattice  sites ,  n  I*  it  is 

continuous  with  continuous  derivatives  everywhere.    A  siE$plc  .th 

these  properties  is 


■M 


•50* 


E  o     J 


is  the  position  of  the  lattice  site  j,  oral  I   li  the  poai- 


tlon  of  particie 


ye  can  dolMuo  a  quantity    Y<(rt.),    •  \M  a  coaposttat  of    iKij) 


a 


*(%)*¥?«<£)  (19) 

j  j  i 

txiat  rsnedber  that  Y*(£)*  aas  aeeai:,    y  nee*     iM  aitoo,  uhidi 

tbftt  2ar£je  values  of  5**(n.  -  v„<)    (jive  r,^)  *  0.   Hence,  ( 

2    l 
presents  no  complication  froa  overlap,  since  ve  have  5  >  >  *~?r  ,  sad 

_^ 
tke  only  tern  Of  r  that  contributes  near  9L    is  r 

tte  redefine  teaas  in  (IB): 

She  subscript  1  appears  on  T^Cl)  since  (20)  on2y  describes  the  density 
near  t**  In  (30),  u-,  is  the  <2in         -  x\J  fi 

r.  to  the  lattice  site  3,  and 

W r> '  ^P  "CiV1**  •  (£1) 

Puttiac  the  aosuaed  solution  (19)  into  (21),  vg  have 


Hating  that  vith    d'*  >  >  -^  »    oil  the  tlowanti  or  the  snsa  en  i  are 

very  sa&ll  impfe  for    t  «•  J,    we  go*- 

*«(\«)  -  (eooot)  /  •  M*j  > 

where    g, *  H,  •  r4f    a  coordinate  with  OBd^n  at  the  lattice  cite 

When  we  write  (22)  la  Ito  diagonal  fom,  where  we  conei&e:*  OuSgf 

rA$)>    w*  fiico  olytaln  only  the  port  of    y{2.)    which  deeerlbea  Its  prob* 

lit?  dene:Uy  iftNMl    IL     and  we  nuot  obtain  the  rest  of    r(l)    aeixsjy 

•o*    I      \)2 
•  torus  with    e  t  f  1. 

Fomuln  (C-3)  of  Appt.  0  telle  uo  that  (22)  ic 


m  I  -oG(u  •  rn)"\ 


r  a 


la  Aj^pondi::  C  we  show  that  any  function    *      in  the  ),  \fc 

ic  aiow^y  teaching  near    iy. ,  »  ti  .     compared  to    c        J      near    q  ■ 
■will  give  ■  see  fee  fora  (23).     lath  (23)  vc  con  Co  ■  prodx. 

of  correlation  lone  in  (20).     The  result  of  the  pfOflUUt  if  displayed 

in  formulas  (30)  and  (31),  which  we  m  oed  to  derive* 

10s    f.     li  well  "behaved  for    tt>  pnft    f-Cu^.) 

about  the  value    tt  L«     (later,    R     ic  to  be  a  distance  to  a 

of  atoras  equ.  B)  arbi  t  paint  co  0  .  •    She 

egpenstrai  U 
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We  define  the  product  («vj  !2  ^t  ,    so  teat 


•i  j,  *•  -~  '~ 


I  e  NM1 
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L  t 


J  J 


*1  f  B 


*     #     * 


We  introduce  0  (juarvtity    a     «herc 


1 


Q     is  defined  ee  follows 
G,   »  in  s 


o  i        A.  a'i 

j    *    ■*■     £.  ■      iff  1      ■  1 1,     <*. 


E  « 


A 


ii  vi 


* . 


(36) 


But    j8n(l  +  *)  •  jc  *  x*/fc  •»■  3tV3  -  .t  A  v  .   ,   .  ,  00  that    GL. 


/^ 


g^  -  j    Aj(f  7f)  *  i&f/W/t)  -  (ApteXf  *)7(?) 


0-  .  fi 


But    (d/luXf'/f)  «fcr  r"  -  ifD/r  ,    00  ve  on  write 
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Ue  partially  Ic^ooed  the  constraint  on    r^    In  evaluating  (83) « 
ife  now  coqplete  the  process  lay  evaluatiafi    GL     for    r^  »  tf,  ,    tdbte  pc 
tion  of  particle  1,  resatiberlng  that    f  *'\*)    is  a  frEw**-**1  of  **» 

distance  frcci  lattice  site    &•    f2tm  evaluation  of    GL.    proceeds  by  per- 

2 
SonalnG  the  Bin    £  &,    and    £  A  ,    for  constant    R.  ,    corresponding  to 
j    J  j      J  * 

a  shell  of  atctas  at  constant  radius*  and  adding  over  shells. 

*,    z  A,  -  l  (luji  -  |«y|  >  -  £  <  ^  *^  -  S%  a^  coal  -  By). 

[Figure  (4)  ohovo  tbe  relationship  between    u,  q    and    ft»]    n;^»adiag,  vo 

gat    £  A   »  £  (Kjjl  ♦  <£^*4lj)  *  %  cos  9^lj  *  **  sA&f^)  * 

...  }  -  WmA*    2he  sua  over    cos  9    yields  precisely  aero  in  a  Dravais 

lattice,  and  in  a  hexagonal  close-packed  lattice  is  aero  in  the  a^ierically 

ayonotrlc  approximation.    Mo  ta&e  tbe  average  value  os    9    to  be  1/3, 

getting 

S        •  S    7^-  (1  -  cos2  9)  ♦  .   .   .  •  8     ™~  ♦  0(<a  ^). 
Puttlag  all  oaabers  of  the  sun  vita  the  sane    i. .    tooether,  \xs  have 


^ 


E  a  -  •:  -4~A-  ,  as) 

-where  there  are    n.     aeabera  in  the  shell  at  dinta.  [see  Fig*  (k)} 

Eliailarly, 

J  *  E  ^^5  +  Q^2  -  20^  %,  coo  9  + 

-  2s£.(i  +  —r-  -  •£=-*  cos  9  -  -4-  cos   e  ♦  .  .  .)) 
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Figure  i+.   Lattice  Sun  for  Two -Dimensional  Square  Lattice 
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„2  coo2  8  ♦  0(c*) 


O(^) 


2 


J  ^ 


S  ^A-  *  o(^) 


tie  did  not  Iieve  to  aosirae  a  Bravalo  lattice  or  spherical  OjOaetxy  to 
elininate  the  tern  MS  0    leading  to  the  result  (2.    . 

We  nov  have  f*as  (27) 


°1  "         3 


1,14 

5*1  B 


.rtf 
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Mb  i*fSi» 


i  .  i  a 
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O.v- 


f 


J\ 


(30) 


Bow,  using  (25)  and  (26), 


rx(D 


♦*U)  ^(i)  p1 

-£t?(l)i(i)  0^«X 

•  E  **(!)  ^(1)     JT  tiPx)^\  e 


(3D 


2hc  tsanipulations  of  foraulttB  (2*0  to  (31)  have  given  us  on  I  - 
preeeioa  for  the  probability  dene  i  i  atom,  deterc&ned  by  correlations 

with  all  the  other  atono  of  t&e  solid*    9di  ixsportant  result  la  Use  pro- 
^ectloa  in  ordinary  threes  ftl  space  of  the  probability  density  de- 

aeribed  in  the  space  with    3W  *  21(8  -  !)/&    variables. 


■ 


011  put  (l£)  in  &  torn,  involving    &,    the  parapet; .-         the 


A     ;;. 


Cauooian  representation  for    r-    Ifc  have    r,(l)  "  «**      *Hir/tf'y/i',    ^r 

hypothetic.    ISeaee, 


2    "y%%S 

i    * 

Store,    A    io  a  nortaalissing  £aetoi*. 

Slaking  the  logarlthra,  ;?e  have 

ia  A  t  |^(l)  |2 
62  -  -  R ~  Sv — -i — *  (33) 

*  3  ^    •     v 

Our  aDOix^pfcica  that    r,(l)  •  o  evidently  Ocpcnde  on  the 

relation 

lia     A  £♦«**(£)  -  c         x    .  M 

^      ^       1    *      x 

'Zto  work  vlth  (3*0>  ««  tssosat  cc,  tlae    $'c>    ia  BUM  detail. 

aaouac  that  they  fit  the  ayesaetry  of  the  aattiee.    For  exosple,  in  a 
oiraple  cubic  lattice,  they  wuM  have  the  Sosa 


coa    O   ,         Uoo    24* 


V     "i        IV      ■        Ua     * 


B  B 

coo  I  t 


/ 


vfcera    n  *  B  *  X,  S,  D  -  1,  .   .   .  6,  >,  ^,  3,  2,  1,  l/B,  l/3,  lA*  •  *  • 

l/H  -  1,  l/ll,  1/tf  ♦  1#  «  #  «     and  the  I  phase  factors  \&ieh  plr: 

n»! 

the  vnveo  in  the  lattice,  ao  tlie-  eo  are  at  the  ssadjaa  for  c^ 

a  end  l/n  and  oo  that  the  nasdUsa  for  oven  a  or  l/n  are  EsLdv&y  between 
atono.  2huo»  ve  have 


•5r" 


a  •  1/3 


a  »  Xfk 


a  *  l 


a  « 


a  «  3 


Evidently  the  vavco  vith    a  <  X    con  be  uoed  to  coaatruct  a  packet  olicrv?- 
iag  the  probability  density  near  individual  lattice  aiteo.    ?he  lone  vavoo, 
with    a  >  1,    do  act  eater  into  audi  pool-eta,  but  appear  %*wa  vra  coaalder 
aeveraL  particles  together.    Tor  exoaple,  ve  caa  see  £roci  the  diagram 
that  If  x/e  coaaidor  h  atcoc  la  a  row,  we  ouat  0©al  vit&    a  »  1,  a  » 
a  »  3,    and    a  *  k,     longer  \awe2ea0tho  vili  not  m 

2a  performing  the  iatecrel  \Q2),  ve  could  £ust  as  veil  integrate 
over    k    aad  evaluate  the  correlation  vith    J,    ©a  perfora  the  integral 
over    di*,#    Shia  consideration,  together  v.-."  faNfetatt  of  the  prodv 

v  f  «(%«)*    leadc  us  to  the  following  interprets  r  /  li*)*. 

It  foxBta  a  packet  of  vavea,  vith  wavelengths  ahoruer  than  about  1/5'"  , 
.near  &•«  Sbereftore,  the  renainiac  $.. ,  which  appear  explicitly,  are 
Oetenaiaed  by  other  considerations  than  correlatioas.    For  eciacaple,  they 


*  \     /  •  \     J  «        /  •  \       % 


■ 


# 


*  •      ^>^      » 


Id  be  caHed  upon  to  penult  ejonstructinc  states  is  to  t 

FkuII  principle,  or  to  xaafco  up  eieeaotates  of  parity. 

Lcudin,1  in  the  third  pe©er  of  hi©  series  on  aany~parfeXe  cyu- 
te»s>  describes  the  use  of  two  such  states*  those  for    a  ®  1    r. 
n  *  1/2,    as«  s^oaetriaed  trial  function  for  conaiderins  pr- 
tions  in  &  system  ouch  as  MMfc 

In  Appendix  H,  -we  emluate  the  sua    S  ?,  (r.)    ia  the  scoeuhat 

odd  nuatber  of  ataas  alone  each  edge,  w^frrirf,i)g  samtaetry  shov      &  .    In 

thio  case,    rl     can  oceino  only  odd  values  end  all  the    8      arc  zero. 

n 

^So  asnuae  that  the  riauli  principle  pemtts  two  jjertiolcs  in  eaah  of  the 

It/i    states    n,  *  1,  3,  .  .  .  (5)^  .     wo  tci;e  the  long  vavdeugth  can- 


ponent  of    ^    to  b© 


♦i<*)  -  ^ HBM&lH B)  •  ;>  - «*»  -  !H)i/3 


\_ 


*4 


Ite  find  that    A  *  l/(  I    c£  k),    and  that  all    H    of  the    C    ,^    ore  etjual 


. 


to  the  acne  quantity    CQ  »  t&?*J  '*<*))}.     In  this  oppression 

for    CQ,    0(a)    in  the  value  of  M  |  .on  product  -,.),    where 

Sj^    Is  the  dlatanco  Ami  lattice  site  1  to  lattice  cite    k*    ^h&  values 
of  the  C*s  do  not  enter  into  (3**),  however,  and  the    n„.     in  the  dencraina* 
tore  of  the  axsuwoto  of  the  cosines  £ivo  the  convenient  result; 

mm 


#r<t  m  c  ar*  «   . 


mn 


Saamaoo  m 


Jvm  have,  for    o*% 
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♦  a    , 

ft2 

«  -  s 
X 

3      * 

(36) 

uhere    ct    la  &  turn  of  orte    if  '*,    hence  M|         i  ui  iter  Xaroo  oyotoas. 

la  £jppea&l}:  X>,  we  derive  esprc.  Of    &.  •  (^ 

8 
bos  two  ports.    She  first  pert,    &Q  *  jepencleat  of  5  ,    i 

8 
another  part,    £vj*    5«!pttwSs  on    &  ,    but  Is  a  aenll  correction  to    { 

la  the  setae  opg?ea&bt,  wo  show  two  e:jpre.  i  for    fcg*    Sic  Bkfjpl 

U 


3  ck 
6*0,1  *  *  ™5 


1  * 


Bo  o 

■  ■■I'd— »y  f»        ».»M«pft 


(37) 


1*. 

vfcorc  B  ^  — ^...£  .    b  »  2/5.60  for  Mr«    Uo  oftteia  (37)  by  apply- 

ine  (30)  to    f'/f  •  (#~)'/^%    wbere    X    is  in  the  series  repsresentatioa. 


S!ie  second  fom  is 


%> 


# 


1  -  x( 


(33) 


-  r(r~r0) 
where  we  use  the  representation    X  *  1  -  (r0/r)e  ,    r  *  0.696, 

s^Eoetry.    Us  have  sho«a  (37)  aaa  OB)  in  ftgm  5. 

!Zhe  pert  of    &     which  <3epen£a  011    5    is 


itr 
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g,  from  series  representa- 
tion of 
g,  -  from  Yukawa  representa- 
tion  of  X 


Figure  5.   Elements  of  the  Logarithmic  Sum,  g^n,  Plotted  Against  Radial 
Distance. 
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[x  *$<&[:;_  -  r0J)] 
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unere    $    is  the  error  integral*    Shore  are  additional  od  is  term 

of  order    ($3^)%    vhich  are  discussed  -in  <ft$pento  D. 

Froci  (35)>  (375  and  (%?),  ws  oan  calculate    6    *»  5^ 

V»  first  obtain    &/"    by  evaluatiae  -^r—  »    a  quantity  independent 

of    &J".    Ws  then  evaluate    t  -~^—    until  a  consistent  value  of 

obtained.    Figure  (6)  shows  the  reoui;,  of  such  calculations*    3ha 
calculated  values  fall  on  a  0000th  curve  x&lch  can  be  fitted,  vith  no 
discernible  deviation,  to  the  fomila 

002  •  B»&  *  fc.SS(a  -  *r.064)2  Jf  (Ul) 

StoaauXB  (&1)  is  cjharacrboriatic  of  Be3  in  the  body-centered  cubic  iatt 
and  is  valid  in  the  rang©    3»*#  <  a  <  3»Qfi.     Appendix  8  Aft  •  seqple  cal- 
culation for  a  point  on  Tie 

She  ■trfgK'paTHi'li  repreoontr  of  ore 

r^(r^5  *  (coast)  E  |^  r^l*"  e 

•  (ocwrt)  0*°  *    e    ° 

-  (const)  e    0<l    [1  *  0{jf  l/3)l 

-  (const)  •-**  •  **■  -  ***>*]  ?  m 
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tm  (J*2)  we  tdao  have  the  complete  eiaele^particle  density 

~C0.™>  ♦  M0(«  -  k.Q6kf}{5,  -  r 
r(l)  tt  8  (coast)  o  J        x      . 


81aple»?artlc34>  SoerEy  Density 
In  thie  eectloa  ue  uee  the  oioglo  particle  p.  /  <3onG 

w©  have  obtained  Itt  the  luat  Motion  to  deteroiao  the  average  energy. 

Ebe  espreooion  contains  the  la.  :jnctant    &    throu$j  tho  probability 

density  paraoeter    0. 

tfeing  operator©  introduced  in  IJq.  (32),  together  vlth  the  density 

natrix  of  (13)>  we  can  write  the    3    particle  onercy  density  aatrlx  in 

tJ3e  Pluvinage  representation: 


<  I  >  -i  /  B  «.  j  jgrU  %(*>  *  £  l?u  [JL 


♦  i  V    (r    ) 

i,r  iy 


1JJ 

-  <i« ;  *»'  iWP  «fieip  -  «#%>  iA»  ug  -  %fti 

d^ai,...  d^.  (U3) 

She  normalisation  integral  is  represented  by    A,    which  3ms  the  fora  or 

the  integral  shown  in  {h$),  with  the  energy  operator  replaced  by  unity. 
She  prine  coordinates  perait  the  wave  furatiUm  to  act  as  a  aon-dla©onai 
density  xaafcrls  under  differential  operators.    *to  aust  bear  in  saind  the 
relationship  between  the  %  and    ..  .are    r^ ,    is  e  dependent  variable 


..•  .  .. , 


useful  ia  &tSmmx%i*tem,  as  explained,  til  Appendix  A.    3be  peaaufcotion 
operator    P*    permute©  the  priae  coordinates  according  to  Penal  ffcHfcUrtti 
<3ust  00    ?    perautes  the  unprinsd  coordinates. 

*Je  sij^piify  (&3)  by  noting  that  the  operator    T^ ,  +  1/2  V       sashes 
no  attribution,  since  [ft  .  *  i/2  V\  J  fcj  »  0    for  ell  i,  J.    5tar  our 
sysxwtry  ve  asaune  that  the  torn  in    n'    is  negligible  coopered  to  the 
torn  in    9?«    As  a:piaine&  mrc  fully  ia  ^peadia  P,  firstly    K'    io  start 
range  because  oC  cfce  telvntive    Mfc&i    "^fcich  vanishes  ct  iar@e    r^., 
sad  therefore    H*    is  proportional  at  nost  to    II    end  not    11(12  -  1)/ 
Secondly,  in  a  Oravais  lattice,  a  3areo  pert  of    XT    cancelo  out  because 
H'    is  linear  in  relative  uUsplaceaeat,  and  a  Given  eton  has  a  positive 
displacement  vita  respect  to  one  aton,  and  an  equal  and  opposite  displace- 
aent  vlth  respect  to  soother  aton»    A  contribution  to    E*    remiss  \?hich 
is  proportional  to  the  overage  OQuare  oC  the  deviation  fron  the  lattice 
oite.    In  Appendix  F  m  present  a  piau;  9  arcuwnt  that  this  con- 

tribution io  considerably  less  than  <^  - />• 


atins  v4(^)   toiiiMtM,»M]««tilait( 


<    > 


sqfe 

la  (¥*)  -we  shov  the  eyctaottlc  perantation  of  all  the  coordinates  in  the 
$*e    in  t2ie  diaGooal  ..trjrcocutation.    ferautationc  in  the  dleoonal  repre- 


floatation  are  the  onOy  oaeo  tbat  affect  the  oinc^particle  operator 
?  .      (Tee  the  iaet  relationship  in  (5)*) 

life  note  that  the  part  in  cur3y  braclseto  lo  the  function 

ex 

in  the  bub  ewer    i    (a  reoult  of  the  oycoetrio  pemutatioa  generated 


S  $«^(i))>  ve  have 
t    t 

\*o  redefine    &,(*w)    *»  its  acra-diagaaal  fo: 
ihat  (hrS  bocoaoo 

Waen  u*  introduce    IT  f.  (u..  ),    ue  evaluate  <^  S  &  >  , 

I         I 

<^  •  ;'.  y  *Z^Z%  y    ,    flince    IT  flLfaO    cteecribee  the  prolmbii 
^      i    *    4 

daaoity  near  lattice  site  1.    2bia  a  an  rooulto  froa  the  oaaR 

CQnoi&eratiaae  on  affected    r(l)    in  Bo..  (10)*    2io  nornalisation  factor 
li 

a  -  /  s  li(D|2  it  &(«•*)£. .  <*3) 

i  &       *** 

Ctaffsytof  Qtrt  tfba  operatione  in  (W,  %*e  nave 


<«%>*rt)!|'^!Vl,W,,j 


i  i 


-  r    '  T^  I 


*  /  «L  C  |t4(X)r  2      rt(t)  :t(r    )  V*  *<rlt)dr; 

3h*  tttocgceto  la  the  second  tarn  refer  to  tlje  non-diagonal  natu 

H  f  .    2nie  tern  of  (¥j)  vanishes  because    V  (r  f.J     contains  a  eici  c 

k 

the  lattice  point*    kj    dotted  Into  the  constant  vector    ^*>  $*  7  $,*     Vpan 

svKsaatian,  cosine  tern  In  tfce  scalar  product  Gives  aero. 

2a  the  last  section  we  learsed  wo  couid  treat    S  |  $«(l)|*     M 


l    hence  the  fixvt'tem  is  e:rpected  to  !»  propo:  to 


^"^  ,    and  therefore  lljgltjitlilii*     In  ^pcattla  B  ve  consider  t  jo 

in  detail  end  verify  that  It  is  proportional  to    ST  ™«     2h  the  flec- 

tion we  eonai&o-*  •  abNogHj  and  find  that  the  o  of  the    $'- 

the  energy  Is  indeed  snail,  though  not  proportional  to    if       .    E  i  or 

-1/6 
II    '       arises  as  a  direct  consequence  Of  neglecting  spin  one-. 

in  (J0). 

She  third  tera  vill  psrovido  the  cohesive  energy.    Hbwcve.  Oo 

the  integral  over    dr,.  ,    we  awot  apply  the  coasts*  p±i       *n  " 


to  the    y.    fuacticsao.    Shis  -wo  do  heforo  integration.    Iho  operator 
will  operate  on  the  reeuit  oi*  the  intecratian.     Here    V.      take*  Oeriva* 

tlveo  of    X(0    at  point  1. 

Sfc  deal  with  the  iafceeral  ^  /  t%M  \t(*lt)  t^it^t  • 

I — 
vc  identity     -wr    X  7^  x    vitfa    F(r,  J,    and  uae  the  result  of  Appendix 

C.    \to  show  there  that  any  iateerel  or  the  ?om    /  rA$)  F(i\,  Jdt- 

given  appro:clrjateiy  I37    ^^i^3-  *  0    c     ~ou         ]    provided  that    F(X) 

is  slowly  varying  near    X  «  u. ,    ccc$>ared  to  the  variation  in    r(3)    near 

2hen>  Oaflalac  the  laet  tens  of  (1*9)  as   <  ?  >       , 

1,0 


<7>i, 
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Bevriting  (50)  and  neglecting  a  derivat 


^*il  ♦  0 
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0Ult 
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1*0 


1 


o  \  ~\ 


5  u 


-T*- 


/*■    WiHM    111,    ■  <!■■»■»  —  *» 
T  *V  ^It^     ~a  alzXttU'v  l2W«tt  to  us  froa  Dg.. 


-  e^*  o-8  *  (const) 


OCT 


n  **(^it)  ■  (const)  ©  ^  .  ja) 

*&  nuet  consider  ooot  cmrefuliy  -the  vegr  in  tixich  u©  ccdbine  the 
contributions  of  tho  Inplaeian  operatiae  •»  each  of  the  f-#  *tei  we 

eoploy  the  xtgweeentatloa  c,  ^  .  Uiia  special  attention  **,aa  also  needed 

•-> 
ia  evaluating  £  &*  «o&  5  A,*"  ia  the  preceding  section. 

Bach  of  the  it*    has  a  siaeie  variable  u~.  or  A,  for 
to  cot  upon.  Each  or  the  texoa  In  the  integral  In  (50)  Introduces 
"kinetic  energy"  of  relative  option  obtained  by  to.  .         (motion 

OOR3&1  to  the  line  between  1  and  ,3  Gives  no  contribution),  and  hence 

2. 

represents  one  decree  of  freedom.    But    v      not  only  tokos  do    .       rea  in 

two  directions  no;33&l  to    EL . ,    or  vould  tract  each  acriber    j    as  three 

Oogreeo  t        -  adon.    Hence,  to  use  tne  op;;  aynpe. 

e"      ^      to  represent  the  influence  of  corrclationo  on  tao  energy 

ataa  1,  ia  the  single  particle  oppsroxtantioa,  \m  nust  toko  1/3  Uie  -tssult 

obtained  by  use  of  the  lapiacian: 


♦  /  C  ♦*(!)  tyl)  »*(*  %2  &0"/2)  y  «*<-  %*  %       )     i 


A-Ze"0"    ■  ;^(i)|2arr  (54) 

it 


In  discussing  Eq»  (%),  ^c  have  already  oeen  that  the  first  tera 
of  ( 53)  is  of  order    fir    and  hence  of  order    13~  '     »    eM  negligible  if 
%*o  neglect  overlap  and  osodbaoge. 

Therefore,   vise  oii^le-partiele  energy,  froa  the  second  tera  of 
(53),  Is 

/  -  >    .   1L    1  0      »   *      sj2  .  (53) 


cchca/je  Energy 

(a)  la  th~  a,  using  the  Pluvlnage  foro  of  the  vave  func- 
tion, we  -be  oyiaaetry  effects  on  the  magnetic  properties  of  the 
solid. 

Entirely  analogous  with  the  treatment  of  electron  e:y:bon©e,  we 
assume  that  the  Pauli  principle  applies,  so  that  ve  are  Halted  to  a 
state  description  of  any  pair  of  particles  <anich  it  ontiessaaetric  in  the 
interchange  of  all  the  coordinates  of  the  pair.  Vq  atteapt  to  find  the 
vave  function  of  two-particle  states  of  ioveot  energy,  one  flya*     in 
the  apace  ecK>rdlnates,  corresponding  to  the  triplet  opin  state  and  the 


nJQm 


other  antleynaetricj  corraeponding  to  the  singlet  spin  state.    9H 
of  these  otatoo  are  ccqpared,  end  we  olrtain  the  e::chanee  energy. 

ETl    is  the  MtagSi        in  state  energy*    E-;l    is  the  triplet  energy.    Posi- 
tive   J    corresponds  to  peraUol  sp&na. 

Our  calculations  vlil  obtain  <^  E  >  (>}  «*  +1,  *l)    and  a  ooamlisu- 

tian  integral  for  each  case,  froa  %&lch  ve  can  construct 

<B>n    «  ^B>  ♦  ia  € 

ftP  *•  ■  W    IWp^wJw       SB         .Li»iW».llill».  .lilt     1%     »i       I 

"n      <i>        <i>  *  j  a' 


Ibon  ve  have; 


-  A^#£) 


J- ^ir~    •  (56) 

Exchange  in  the  may-body  caoc  is  souevhat  difJflerant  than  in  the 
two-particle  case,  for  which  we  have  Just  defined    J.    "We  define  e::chan®e 
as  exchange  with  respect  to  a  class  of  neighbors;  for  erxcsple,  exchange 
energy  with  aeert  nearoet  neigibors. 

He  will  conipare  tfce  energy  of  the  solid  in  two  cases— one,  when 
the  spins  of  1  and  the  class  of  neighbors  are  parallel,  and,  two,  when 
these  spins  are  antlparallel.     If  the  e:;ehau£o  energy  with  raopeet  to  a 
given  sat  of  nelGhboro  is  positive,,  that  e::chaa£e  integral  is  teraed 
ferrceaegnetlc. 

Vto  proceed  by  writing  down  the  Itailtoniari.    lie  delete  the  tara 
[X/2  V,  4  +  t  J  because  it  yields  sero  In  the  product  with    X,  ,«    Accord- 


t£g»o»  «C     .w.r. 


HI  M*»J  .   —W      V:o 
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ing  to  Appendl;;  f$    li\,    will  appear  oeain  as  a  aymotxically  oriented 

eua,  with  the  linear  velocity  dependence  giving  aero  to  lovect  order,    tfc 

ttifrft    V^Cr.)  to  be  a  constant,     3m  Homiltonian  (for  lowest  flrtUW  ngi 

-^y)  Is  then 

H»ST.    . 
1     1 


We  arran&e  ttWN  operators  Mi  follows: 


1,A 


(57) 


Here,    X    includes  oil  oC  a  class  of  iMicnboro  say,  ms-rt 

neareat  aeiehlors  of  the  acne  spin,  vhieh  cannot  bo  distinguished  vithla 

iQ/naetxy  of  space  and  spin*    V?e  have    I    K '  1  ■  I*  -  1. 

X 


Ite  nov  set  15)  the  eiiergy-clcaaity  aofcrix,  ao  la  the  prcced 


$0      '■    &    •    :. 


1  M      ^ 


li,x 


A     (1*)  d)     (2')  6     (3*)   -  -S 
^  (1»)  *b  (£')  *b  (3*)  - 
1')  *e  C«f)  *c  (3*)  - 


saae,  unprla© 


±  I A  VV  w  %  •  *  •  %  • 


(53) 


is  a  mi—1 1  mil  1  Mil  tlitagml*     note  that  a  ehazice  in     .    only  chances 


4.      VJc  itereated  in  setting  up  ftjMNftf£Mi  states  betuoe::  1  S 


IB  I 


A  r.v. 


J  flfif 


all  the  other  particles.    Ml  have  no  interest  in  the  cyaraotry  of,  oay, 
id  3  as  a  pair,  since  vc  Co  not  propose  to  disturb  their  cgraaet.    i 
Ue  use  a  fcm  for  $.•     bore  has  a  natural  ayaaetry  for  I 


(t) 


7DC 
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each  wave  has  a  phase  to  give  the  creates*  sycEaetry,  with  no  nodes 
lattice  points.  Equation  (59)  is  written  for  a  ainple  cubic  lattice. 
For  other  lattices,  •     p  functions  can  be  constructed.  ?0r  exaapie, 
in  a  HaCI  typo  lattice,  tfeere  the  Ife  represents  one  epin  aystea  and  Gi 
the  other,  vc  would  have  two  oeto  of  equation  (!£>)>  «ech  with  U/ii  aec*- 
bers,  end  each  \rith  a  ~>  (2/r3)a,  the  length  of  the  cube  edge. 

*te  con  rewrite  (9&),  ueinc  a  Slater  deteraiaoat  vhiah  contains 
only  inforpation  about  tfte  poire  1  -  n  in  ontisyaiietilc  states, 
which  contains  no  teras  Giving  ontlsymetry  to  the  voire  3  -  n,  (J  j*  u> 
aj^i).  2h  rowrltinc  (50)  we  oust  include  the  other  ataaa  la      ade- 
quate to  construct  the  contribution  to  the  probability  density  of  the 
solid  fraa  the  pairs  J  -  n.  '    .do  this  by  Including  tbeti  as 
•  !$*(&)!''  (4  J*  !)•  **•  asa  use  nant  foznat,  putting 

■WQ  for  all  the  tenas  not  involving  overlap  with  aton  1.  For  esonple, 
for  throe  particles, 

eU)   b(l)   c(l) 

a(2)  bi     0 
a(3)   0    c(3) 


^(1)  £(8)  4J0) 


X)  4£(2)  0 


♦?« 


£(3> 


■  A*  -  B+  -  C*)(A  -  B  -  C), 


ta£ 
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vfliom  a  -  ^o(i)  ^(e)  <Jc(3) 

3  -  *a(2)  ^(1)  <ie(3) 

If  y»  edopt  the  convention  tlint  vo  drop  terras  shovinc  overlap  of  tfe 
atone,  tfce  product  ■UpltffH  to 

A+A  *  3*B  ♦  C+C  -  A*B  -  A*C  -  B*A  -  C*A 

tpfil^a>l8IVa>l  HfcWI* 

-  (**(3)  *c(3)  [**(D  £(«)  *ft(2)  ^(1) 

*  *n(i2)  *b(X)  +n(l)  V2*3  *  (oB:,»'   ~*3'  3  ^2))* 

Hence,  v©  con  uoe  the  dote ■      P  to  consider  eixhanGe  of  1 
end  2,  or  1  and  |*  but  not  2  and  3>  nor  e^rchange  involving  taree 

Ida  m  . 


.  tl»rofore,  vrlte 


It*. 


Wi,J    W*'J    *a(3')    *a  (*'>-- 
(*•)    ^  (2')        0  0  — 


♦c(3') 


♦*  U'>      o 


\wBm,  prlxie  -*unprxiae|     IT      X±,(:  JOi^ 


' 


►ft(V>    - 


[60) 
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J£  ve  understand  that  vc  include  overlap  between  1  and  only  one  other 
aton  at  a  tise,  for  each  Beaher  of  the  clasa    %    of  tod  .  raishabie 

aec&ero,  (60)  gives  the  ocoe  result.    Using  one  of  the  Mtf 
out  of    %    to  xtjpreaoat  all  in    >.,    ue  | 


5?,  2L 

«i  *  r  n       -ii 


V1,} 


%Tv*l 


V8'} 


0 
0 


^(3')         ♦^•)- 


o 

i 
i 


(3*)      0 


»(V) 


"3 

i      \ 
i         \ 

\  \ 


|  oaae,  prlae    «•  unprlm  |      I!      X^*^)  W'ij^  •  *   •  ^g  • 


'- 


We  calculate  the  average  value  of  the  operator©  using  a  dif.  netlon 

each*  hence  uc  have  a  different  noroalisinG  constant  ftff  or 

lte  see  tfefft    Tx    is  cpi&c  "to  contribute  <  f.  ">   to  the  sun,  the 
.ue  obtained  for  the  direct  energy.    We  can  evaluate  (6l)  far  any  vain* 
of    A    foroaJJy,  and  obtain  a  fornula.     vJe  can  apply  the  taM&ftj 
for  each    X,    adding  the  results  to  ©et  the  total  energy  eboaot  I: 
iac  the  spin  orientation  of  1. 


\*-<«i>-r*;/4>\ 


Ttja  ¥rjn 


cceao,  priae  -*uuprtoed 


CR 


(6s) 
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Iii  ((&)>  v?c  liave  aood  only  that  port  of  the  cJetcrninaut  \jfcieh  C  ooe 

the  oyupotrlaod  vave  function  of  the  pair    (1X)» 

Vfe  r*ow  adopt  the  convention  that  the  Xnfl&z  identifying  the  ixt?ia 
of    a,  tj,    is  plus  one  ft*  the  spin  of    a    opposite  to    A#    and  minus  one 
the  spin  o:    i    parallel  to    A.    2hen,  the  concession 

is  antigyaaatrie  In  interchange  of    1    and    A,    end  the  rtauli  p 

to  otoe^a*  if  ve  use  (635  ia  conjunction  with  any  s^raractric  function  of 

X    and    A. 

Yte  am  rovr,'..:         )/Mipi'0>«fit1ng  all  atone,  except    1    Ml    A, 
1  sin^ie-partiale  density,  as  ue  i  (35)  and  (^5): 

•  *|*A  ' 


She  factor  2  in  {&$)  arises  feftosswo  of  the  c  «t  the 

nance  tOHM/  is  tiso  dlfAua  »  two  particles  upon  exchange, 

ubareas  ((&)  t&ves  the  difference  for  one  par.  SorsaOa  (Oh)  is  the 

result  w  wuld  nave  obtained  If  tjc  hod  oosiraed  ttW         KWB  1  cad  atom 
give  sri  a<chanQe  «nesey    £*gf    t3>ea  ^°  toto^  o;'.chaacje  obc  all 


-fC- 


i»ie»bors  or  type  B  fj    a^  J^>    *«&ere    a^    to  the  nuober  of  type  2  otoao. 
However,  under  uh.-s  sickle  MVttqpt  the  question  arises:    hov  oea  vo 

include  the  I  eaercy  of  atou  1  sir^tamcously  ia  all  the  pairs v 

fQnei  derivation  {(&}  anovero  this  question  by  iKNUf  ■  that  the  klaift 
energy  of  o-'taa  ~  «*oea  not  affect  exchange  at  all. 

She  evaluation  of  (6k)  is  a  aajor  tasi:.    Mo  can  use  beanaaic  oscil- 
lator vavo  function*}  for  the  $*s.    Sfcey  will  include  tho  cut-off  off  act 
of  the  product  over    h    at  lareo  distances  frets  the  pair  (PO,  and  viil 
peroit  the  eymotry  to  be  affected  by    tj,    if  ve  take  one  of  the  ctateo, 
cay    A,    to  be  tho  first  edited  state,  and  the  other,  say  $  }  I  uhe 

NBA  state,  in  a  variable  describing  the  dcGree  of  frecdoa  along  **» 

lino  ^oiniae  #*o  lattice  sites  of  the  pair*    !2uc  function    IT    can  be 

«ir*  r    ** 
reprooeated  as    IT  *  (1  •  A  e         ^  ),  A  »  !.$&,  b  -»  0.0  >* 

handle    $    and  in  the  way  described,  they  yield  elewontary  Gauoeiea 

integrals  (with  oox^lications  introduced  by  the  cut  off  in    >r)»    Intro- 
ducing the  UJhwr1]Hi'  of  the  product  over    k    in  the  region  between  the 
lattice  sites  introduceo  a  crucial  problea,  in  itself  worthy  of  a  neper: 
study. 

Instead,  vo  calculate  an  apprcciinate    J,    va  lov  flgiiiHti«j 

using  the  aln£le-particle  density  NtfHtt  . 

lie  learned  in  the  last  r  i  that  r.tf*  could  gtft  <  8  >  by  putting 

rCl'A)  *  M'V)  7(1)     {provided  only  that       \T.  I  oCe^^'^'Su)    van 
otmil).    HB  shall  use  this  idea  to  reduce  the  de  ..-•'*}  to 

the  single  particle  fom.    Ife  can  interpret  the  wave  function 


TtSJQB 


.«r  uw 
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((Ja(l)  4A(A)]        rT/  r<k)  3?(lk)  ^i0l§£~   to  be  the  solitude 

for  particle  1  to  be  near  lattice  site  a  with  particle  h    near  lattice 

site  IU  Shis  is  because  (a)  the  product  over  k  keeps  the  perticlee 

sear  these  two  lattice  sites,  and  (b)  **(r^)  prevento  them  from  being 

on  the  eaae  lattice  aite.  Ml  recall  that  a  representation  lor  x(l)    was 

-6~(K-r\)2 
Se         .  We  want  to  build  up  a  non-diagonal  tvo-particle  density 

aatrix  from  thio  representation,  and  atill  have  a  result  vhich  shove  the 
neoeosary  conditions: 

r(i2|i2)  *  r(D  r(£).  (66) 

We  recall  that  ve  had 

r(D  -sr-U) 

r.(l)  -  (const)  e         J  . 

2he  antieycwetric  fona  of  the  density  netri:r  which  we  seek  is  i 

r(2)  ■  UrJiTrfe7)    ♦  n    ^^JxTrjJ7)  ) 


a(i)  rb(s)    ♦  n    ??J2)  rb(D 


I    ffr ^^^^^H 
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+  n  e  J 


ij  e  J  itmao,  prine  -*  taiprliae} .  (6?) 

Fonsaila  (67)  yieldo  the  conditions  (66),  and  given  the  ooiif 
single-parte.  go.    3*»»ula  (6?)  vould  over-estirate  *Ae  I 

ve  were  to  use    -  h'/v&i;  for  the  energy  operator,  in  the  asae  -way 

as  la  Eq.  (5l)«    Ws  easst  introduce  the  sene  tfector  of  1/3  hero. 

shall  use  the  foxn  in  (6?)  as  a  first  appaxxKixaatiori  'mild  only 
so  Ion©  as  the  spherically  syuaetric  approjiimtion  has  little  effect  on 
overlap  of  adjacent    r's«     2be  spherical Jy  oycsaetrle  ftpprojdnai 
poor  at  saall  etceaie  separation©,  so  our  results  can  be  expected  to  Jr. 
validity  only  at  lower  densities. 

We  can  rewrite  (67)  as    I*  '  «  |f  |  j 

/  -  v  it*-!-*)1-     -^)23     -  v  u%w  •  •     ii 

e     ■      ~         X  *      t  (63) 

Here    b  «  a/£.    She    s    axis  Joins  the  two  lattice  sites    1  and  2,    end 
a.     and    I  ■  nj&asured  from  the  nidpoint  of  the  line  Joining  the  ll 

tice  sites  (that  ie,  there  are  lattice  sites  at    s  «•  !  b  a  I  a/. 

Xdentl:ty±ng    2    with    A    in  (6fc),  and  writing  the  second- 
density  matrix  there  la  the  fora  (64),  ve  g 

En,2  '  /  *#  (-  s)      f  *  tei  *a 

and 


tto  get,  perfona&ng  the  operations: 


r™>  Jl  -?•  ^  »2(2) 


tf~     i    a4*  -  n[2tr  y  -  s'*3   e 

fO        Wt»- **■*■  p  Mm  i      n     HI  imii^|ifciw»>,       III  ill  I   I  —It  UniMWPMfia— ■■iniii-i^<w 

£5  I  ^rrr^ 


Hence,  fi»  Ife  (5$),  vo  have 

-  list     t  ^  Jt    t.2,     -  a   d  a 

■  -    3,2  •  SU  l^    o    -  &  J    e 

a  *  e 

1  -fi 

<  17-  1* 

W»  uoe  fiq.  (56),  but  include  the  factor  2  from  (65) 
x2(fi  -  ^4#^> 


J     a ..,,    ...■„.„. i.T^i     *~~- 

<1>  - 


cting    Atf/<  1  >    coopered  to    <  I>,    ve  bave 


^  -■     t 

J2    tf~  cT 


*  •  -  %  r  *tf-   •     '  (6p) 

Equation  (6j,0  cannot  be  ©injected  to  "be  valid  in  the  hi$b  aeneity 
it  1  since  the  overlap  of  the  functions    f(k)    entero  strongly  ia 


-Bo- 


il bettor  egpp&mi&iBtioti  would  conoicrt  of  uoincj  less  tfcMB  opk*. 
igattftry  in    r1(i)>    fox-  vxma&la, 

-  ftJw  ♦  y"  ♦  **)  *  B 

V-U)  «  e      - 

Ifouever,  if  ve  anticipate  that  these  correction©  to  the  exchange 
integral  are  of  the  omm  oign  aa  (&>)>  tben  (Oj)  gives  a  lower  liiait  on 
J. 


ran  v 

Elastic  Effects 

■■ —       -  -  -  i  mi  a    i 


.-«*..»._     -   ._- 


In  Chapter  IV,  wc  bsftt  obteinc  ■         -  -p-jession  for  the  eofceu. 
energy  of  cc  lie  body-centered  cubic  pt»c 

I  CD 

With 

a0  (•  -  k*  "    . 

3he  ejqpreosioc;  for    &0      is  obtained  by  const  -g  the  ground-cstete  wave 

function  for  the  elven  lattice  structure  from  the  correlation  t 

,c/r),    gftlmiletiaft    5Q      for  various    e,,    and  then  fitting  a 
curve  to  a  quadratic  equation  in 

Foxaula  (l)  gives  the  cohesive  energy  directly. 

i 
Goodfcind  end  Palrbeaa,    end  itelon    have  observed  the  temperature 

dependence  of  the  nuclear  resonance  relaxation  tines.     They  use  the 

•a 

Sl&ertbergen,  11  end  y  of  the  diffusion  r 


i     i      >odfclad  and  V.  M  ,  Helium  OSiree  (1  by  J. 

Daunt,  Chio  "tatc  diversity  Press,  Colura&ua'J  dnio,  S55)  p. 

lited  by  J.  0.  Daunt,  Chio  state 

Univerr.  3,  ISTio,  1,60)  p. 

.^scfcergen,  E.  Purcell,  en  ys«  .*v.  72,, 

3)* 
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anism*     The  relaxation  time  Low  the  Arrhenius  law,    T  *  TQ  e  , 

where    E.     is  the  active:  ergy  of  dlf.  -In  ,  ve  have 

plotted  the  c;:peri$3ental  da*  the  activation  energy 

gcther  with  the  cohesive  energy  glfM  by  our  theory. 

To  a  first  approximation,  the  cohesive  energy  is  the  I  to 

create  a  ve  Once  a  vacancy  exists  at  lov  density,  little  energy 

is  required  far  the  vacancy  to  move  between   iem'ect-neighbor  sites,  be- 
cause the  next  nearest  neighbors  do  not  Intrude  into  the  space  ia  ,uhi 
diffusion  ta^ec  pit:         H»;?ci,  the  diffusion  coefficient  (e  relaxa- 

tion tirseo),  should  be  detc  .,  in  the  ualn,  by  the  number  C  "incies 

ttoomt  the    bltzoenn  re  Btlonetdy    I  ■  :;,  c--\-  \      /{!:', }.    Warn    L 
is  the  energy  to  create  a  \  a  a  dilute  system  of  vacancico.     Then, 

D    «  lj     D  -  DQ  eap(£      /(:  T).     The  4icatee  general  agresner 

between  toe  diffusion  and  our  theoretical  values  of  cohesive  energy, 

if  we  assume  that    F,    *  .At  high  densities  the  diffusion  active- 

s 

.  energies  lie  w  e  our  curve,  indicating  that  our  aechenisa  of 

ahoul  ude  a  density- icpendeat  activation  energy.     Our 

theory  cc.  rcct  this  problem. 

We  have  attempted  to  predict  the  transition  to  a  hessagonal  close- 
packed  structure  at  high  densities,  using  our  spherically  oyzanetr 
particle  At  This  is  done  by  calculating  the  cohesive  energy  in  the 

two  structures,  the  structure  wltb  xowest  energy  density  being  the  more 
stable.     Our  results  are  equally  valid  far  hexagonal  cloee-p©eI;ed  and 
face-centered  cubic  la  the  spherically  qymetrie  approximation.     This  is 
a  better  approximation  for  the  J*8-fold  syesae-  the  face-centered  cu! 
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lattice  (vith  respect  to  nearest  end  ne^  nearest  Md#i 

6-fold  syuaetry  of  the  herjogonal  lattice,  .nd  the  cXceo-paclJed 

lattices  to  be  less  atahlc  than  the  body-centered  stiMe  for  all  &  015 

.y  to  cxpcriaental  f iadiuea  of  a  transition  to  hexagonal  clo: 
poshed  structure  at  high  01  However,  the  spherically  eyasac*'  s- 

paroxixaati  1  is  to  &U  tat*  between  he^a^onal  close-packed  m 

face-centered  cubic  lattices,  and  it  ie  therefore  necessary  to  put  the 
syaroetry  of  the  lattice  into  the  single-particle  density  In  order  to 
zaake  a  theoretical  stisSy  of  the  transition.     Because  of  the  expected 
eicdlarity  of  clooe-paciced  wfcla  and  spherical  ayaraetry,  we  find  the  face- 
centered  cubic  lattice  to  be  excluded  at  all  densities ,  cot^ared  to  body- 

•tofemA  obWj  . 

2fae  pressure  is  given  by 

P^oV  • 

Hence,    p  m  .  ||    g|  ,     But    V  »  a3/l.30    A3/atoa    for  a  body- 
centered  cut        1   sice.    Iherefarc, 

9  -  .  »  (k) 

3*     — 

la  Figure  0  we  have  plotted  pressure,  based  on  (X)*  ;2)  aau 
Also  plotted  ere  lov-tenperaturc  e:cperiiaental  points  reported  by  B3Mards, 
Sana,  Brewer,  Daunt  and  SfcVllUaao    and  fydoriafe,  ly.' 


0.  0.  Edwards,  J,  L*  Baiaa,  I).  P.  Staeuor,  J,  0.  Daunt,  sad 

McWillistao,  Heliua  Hhree  'Edited  by  J.  G.  Daunt,  CS)  ..y 

:;s,  &>ludbus,1££o-,   1>60)  p.   1 

•ydorlefc,  &.  L.  Hills,  end  1,  -     iy,  Phys.   ,sev. 

^  *09  (1^60). 
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Bernardee  and  Prlaafcoff  (3. P.)  do  not  give  a  theoretical  value 
of  the  pressure.  However,  the  tabulation  of  cohesive  energy  shown 
ftable  II  of  their  paper  elves  on  abrupt  change  in  preooure  froo  positive 
to  negative  at  a  lattice  constant  of  about  3.7  A,  vita  a  negative  pres- 
sure at  hi#i  densities.  Sueh  behavior  is  not  consistent  with  the  emll 
aaount  of  o;3>orlaentaI  data. 

Vte  con  calculate  the  Debye  temperature,  0,  throu$i  the  bu. 

7 

compressibility  usine  an  approsiiaatioa  enploved  by  3emarcJao. 

2be  bulk  aoduluo  io  related  to  the  elastic  constants 

3  -  v  g  -   (cu  < -2  e^/3.  (5) 

W»  cannot  determine    a.-    and    *•*    separately;  however*  Mi  con  aoavsae 
that  we  have  isotropy    (cj,  -  fi  cUi^>    ^  tSlftt  tbe  C****^  relation 

bolde    (c^  *  c^).    Ut  than  hove  3  •  (c^/o)^*    itoere    0    is 

tfce  density  and    e-    and    ct    are  the  longitudinal  and  transverse  wove 
aound  velocities.    »U»Im     I  •  9  •„/>.    B*  M*.  t^mtur*    0 
is  defined  by 

8»hc(6  *2  H/v)^3,  (6) 

where    3/c3  •  l/c~  *  £/ct3.    One  therefore  gets    c  *  0.6^  cM  . 

tte  now  cctibine  the  ftwwtulfMj  assvBagptione  of  Bornordes  with  our 
value  of    E    to  calculate    9.    Ut  have 


Tl.  Bernardes  and  E.  PriaaSjoff,  Jttya.  r*v.  lift.,  1^60  (ly6l). 
'11.  3ornardee,  Ryi  120,  GOV  0)« 
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6  *  h(0.64) 


5p 


ul1^3 


V 


But         p  -  a(1.3/a3>    *V*  -  X-3/a3,    *  •  ▼  **§   4   *4   £ 

3a  3&" 


I  3&  / 


2    dE 


Coobialng,  ^e  get 


Knot 


8  •  1* 


.36  vC 


E  »  0*99   *  U.22(«  - 


2-J  ■     (a  In  A), 


(7) 


9  •  4.86  /fc. 22  fa  -  |  (a  -  4.064)]2  •  k 


9m%Zfr  (oin  A). 


(8) 


Ihcre  ore  no  experimental  data  giving  the  specific  beat*     However, 
near  the  celtinc  curve,  ve  predict  a  Oebye  temperature  of  29.2°  IC,  under 
the  assumptions  concerning  the  transverse  modes  which  ve  described  above. 

Ve  can  also  check  our  results  vith  the  findings  of  Dugdale  and 

8 
MncDonald,    that  the  iindeaann  melting  formula  gives  the  same  value    A 

for  all  inert  gasea 


A  • _-£ «  1.14  x  lo-3 


Dugdale  and  D.  K.  C.  MsoX*»gJj&,  Phil,  Mag.  Jj,  611  (1 
(footnote  4,  Chapter  II). 


-  .  - 


M  ■  molar  oass,  groso^aole 

IL  «  molar  energy,  oals/jaole 
9  •  Dsbye  temperature, 

Ha  can  calculate  a  value  of    A    for  very  low  temperatures. 

taise    a  *  3.6  A.    Using  our  value  of    £    and    9    for    a  «  3.3  A,    ve 
obtain 

A  «  £,6  1  ID3. 


Tbis  value  is  about  tvice  that  for  other  rare*&as  solids,  'low- 
ever,  a  factor  of  1,*  In  8  t  could  cause  the  discrepancy.  Further,  HeJ 
has  extreme  quantum  effects,  and  (?)  any  not  hold. 

In  Chapter  IV  we  have  obtained  on  e-xpreesioa  for  the  e:^honge 

•a 
energy  of  solid  lie  la  the  body-centered  cubic  configuration,  expected 

to  be  valid  at  lover  density, 

2  2 
$    Jk  Jt       •   9, ft.., 

J  ■ "  5  H2-  e     •  (10> 

rIhe  er-perinental  data  of  i'leyer,  Adams  and  Pairbank    can  be  used 
to  eatiSKite  the  Curie  temperature    1L    as  a  function  of  lattice  constant. 
Iholr  data  ore  shown  in  Fig.  0.    We  estimate    Tc    by  using  the  ."  -ion 

?c    contained  in  the  Curie-lJeiss  lav 


*B*  ZKrigit  Adcms,  H.  Meyer,  and    .     .    '. stSfeMk*  r&llua.!Ehree 
(Edited  by  J.  G.  Daunt,  Ohio  Htate  Univer:  .-ess,  Columbus,  Ohio,  1S&>) 

p.  5?.    2he  author  is  Indebted  to  then  for  raaklnj-  their  data  available. 

J.  I«  Van  Vied;,  J.  phys.  radium  12,  26c      I      :.). 
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Figure  9 
SUSCEPTIBILITY  OF  SOLID    HE3    FOR  VARIOUS   PRESSURES  VS    l/T 


*  -  ^V  »  Tc  ,  (11) 


viiere  x  io  the  oseaetic  bulk  susceptibility,  C  is  the  Curie  constant 
for  &?,  1  io  the  teogpereture,  and  Tc  Is  the  Curie  teotperature.  'Ihe 
Curie  temperature  io  given  lay  the  exchange  onercy 

•> .   ;i 

.   ft,;„ft  . 

A 

Sere  J*  is  the  e^ehaaee  energy  (ft  aegative  nuraber)  with  reopen t  to  the 
class  of  neighbors    a,    and    a*     is  the  maaber  of  >  labors. 

t-Je  have  oatiisated    Tr    for  various  values  of  a  from  the 
swatai  datn.^    *fe  obtain  the  values  by  rcplottiae  the  susceptibility  as 
l/x    agBi.nct    T.    Ihe  lov-pressure  curves  thus  plotted  have  an  unoiot- 
able  straight  Hm  portion,  uith  a  necntive  intercept  on  the    t    axis. 
Ihia  intercept  gives  the  value  of    Tc,    provided  that  the  tetapavat. 
of  the  portion  of  the  curve  fron  vniah  the  extrapolation  ves  siade  lies 
well  above    5L«     In  £lg.  10  Me  show  such  a  plot  for  c  35 •? 

etcu  la  Fig*  11  ue  shov  the  results  of  several  such  plots  fTCH  the  date 
of  Fit*     ■    rihe  e:rtrapolation  process  introduce  unoartalttty  which  we 

have  indicated  on  the  figure.     C&  the  sane  figure,  w©  have  plottec 
obtained  from  oar  the*.  lie  agreeeaent  with  experinent  ia  good  at  large 

values  of  the  lattice  constant,  but  obviously  breads  down  for  denser 
solids.    Since  our  theoretical  value  of    J    ignores  the 
syncsetry  of  the  locations  of  nearest  and  ne;it  nearest  neighbors,  a 
not  include  the  effect  of  m  areot  neighbors  on  Mbfl&jtj  one  would 

be  led  to  introduce  cubic  synaet;y  Into    y(l)    instead  of  spherical 
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0.10 


TC(*K) 


0.05 


a 


3-55 


* 


Extrapolated  from  Data  of 
Adams,  Meyer  and  Fairl 

a  68.0  atm 

b  6k. 6  atm 

c  57.8  atm 

d  Vf.6  atm 

e  35-7  atm 

Theory 


3-60  3.70 

Lattice  Constant,  a(A) 


Figure  11.   Curie  Temperature  vs  Lattice  Constant 
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ojptaetry.  -^bexo  would  then  be  e  hope  of  predicting  the  sharp  brea:  in 
2L  "below  a  *  3,6  A  as  a  result  of  introducing  interactions  with  ne;:t 
nearest  neighbors.  Calculation  of  J,  using  Eq.  (IV-6*>),  would  also 
piredlct  t..  of  next  nearest  neighbors. 

We  can  construct  a  picture  to  c;nploin  the  experinental  data,  as 
follows: 

(a)  At  large  separations,  correlations  are  related  and  there  is 
coes>arati  vely  large  overlap  of  tbe  atcc&c  wave  functions,  iho  exchange 

cgral  is  characteristically  negative  because  the  lattice  structure 
end  correlations  tend  to  hold  the  atone  apart  in  a  region  \Aa. 
<*V(x\  J/s^,!  is  positive.  2&e  syEsaetric  space  configuration  p. 
the  atcoB  to  he  closer  together,  with  a  lover  potential  energy,  '-he 
kinetic  energy  is  changed  but  lit*. .la  between  tfee  states  of  differi: 
eyaaetry.  (Uote  t2iat  here  we  do  not  rely  on  the  wave  function  %    to 
carry  the  interaction,  but  refer  to  the  interparticle  potential.) 

(b)  As  the  atoos  are  brought  closer  together,  correlations  be- 
came "titter,'  the  region  of  overlap  between  nearest  neighbors  beecoea 
snail,  and  the  exchange  integral  likewise  becomes  onaller. 

(c)  At  still  lower  values  of  lattice  constant,  the  correlations 
between  nearest  neighbors  become  so  "stiff"  as  to  greatly  reduce  ■: 
however,  correlations  between  Mact  nearest  neighbors  are  not  as  wm 
tive,  and  appreciable  overlap  (and  also  exchange  energy)  vili  e.dot 
between  then* 

All  the  preceding  appears  to  be  substantiated  by  the  "behavior  of 
our  predicted  J  with  changing  lattice  constant.  .1  the 


picture  without  vl  n&  *bat  we  have  learned  fraa  the  theory - 

coasaun    *      for  an  oi  ferromagnetic  siihsteaoes,  with  Oe-.-roaolog  tetsperature, 
to  show  a  teopoxetum  dependent  iaex^eaae  of  the  quantity    XT  It 

coaetaj.it  value  of  the  Curie  3aw,    W  +  3*    She  quantity    KB    MMlHMi  a 
mxicaia  at  the  ileal  tesperature    T«»  0NK  tcc??eraturcc,  la 

to  **ro  line:  This  behavior  is  predicted  lay  'the  theory    '*'    vhen 

interactions  with  next  nearest  ■  are  &  ft  a 

structure  where  alternating  adjacent  crystal  plane©  carry  IdtnUoel  spins 

each  piano,  but  opposing  spine  from  plane  to  plane. 

IB 

Furtfce  fcBOWl      that  when  fee  suMatticee  retain  their 

^relations  (e*g.#  nearest  neighbc  -lei,  neat  noaiT: 

;hbore  parallel)  that    X(f  *  0)  »  2/3  x{T  »  1L).     In  this  express! 
1L    is  the  Seel  teapemture.     If  there  la  ■  "between  sub!* 

tices  (o.c.,  nearest  neighbors  without  correlation,  next  nearest  neighbors 
antiparallel,  or  ordering  of  che  second  Iliad) j  this  relai         Hip  does  not 
hold. 


.  J,  Gorter,  Proceedings  off  the ^Sgtfflpaiva  c: 
'iree  (Chio  State  bniverslty  nn: 
tosearoh,  .'.;?:  p.  199 « 

"She  Heel  temperature    T„    la  the  tesaperature  beicv  which  € 
entif anemone  t  i  shows  mgnf  -darin*,, 

netic  field.     She  Curie  teqper&ture    Z,    is  defined  as  the  parencter 

in  the  Curie-Weiss  law,  Bo.,  (11).     In  general,     fL  <  T  ,    with  their 

(never  greater  than  3)  detenained  by   .'■      elative  magnitudes  of  the  t 

.teraction  between  nearest  neighbor©  compared  to  that  between  m 
nearest  aelgibore.    Van  Vleck10  and  T.  ffegBsaiyr.,  K.  Yosida,  1 
Advances  in  i  k,  1  (1^5)  (:■■  cm«- 


ite:Terence  to  9I(    ;hows  that  the  data  for  68.0  Atm  have  a 
definite  maximum  at  X  »  (const)  x  8,6,  T  »  0.1°  :;.  Ihe  t^uMfl  value 
of  X(T  -  0}  if  3L  •  0.1*  I  is  then  X  »  5.73  (const).  The  observed 
value  at  0.0625°  K  is  X  *  (const)  5-B>  end  falling  slowly.  Hence,  the 
data  at  60.0  Atm  night  correspond  to  an  antiferronagaet  with  a  Her 
temperature  of  0.1*  K,  with  interactions  between  nearest  neighbors  and 
next  nearest  neighbors  of  equal  iegportance,  conforming  to  the  theoretical 
prediction20  that  X<T  •  0)  -  2/3  X(T  *  ?a). 

At  higher  pressures  the  overlap  between  nearest  neighbors  will  be 
reduced  still  further,  and  interactions  with  next  nearest  neighbors 
should  predominate.  *2hea  the  relationship  X(T  *  0)  •  2/3  X(T  •  T„)  will 
no  longer  hold.  !Ihe  data  show  no  tendency  to  follow  this  relationship 
at  high  pressures,  supporting  the  conjecture  that  the  exchange  interaction 
with  one  of  the  classes  of  neighbors  is  very  weak  compared  to  that  with 
the  other  class  of  neighbors.  Our  exchange  energy  becomes  small  at  high 
density  because  of  small  overlap,  resulting  from  the  strong  exclusion  of 
the  correlation  function  XT  at  small  distances.  It  can  be  e:qpected 
therefore  to  represent  interactions  between  nearest  neighbors.       it 
continues  to  fall  as  the  density  increases,  we  may  state  that  the  theory 
indicates  that  J  (nearest  neighbors)  is  small  at  high  density.  Slien  the 
observed  marked  increase  in  Tc  at  high  density  should  stem  from  an  1 
change  Interaction  with  owl  nearest  neighbors.  A  magnetic  transition 
might  occur  when  the  ne:rt  nearest  er^hange  integral  grows  largo.  A  spec, 
heat  peal;  at  the  transition  temperature,  which  might  be  a  ccsaparatively 
high  temperature  at  elevated  pressure,  could  explain  the  difficulty  in 


.  O ,      •  -BO 


iom**<  I 
ttff 

ft* 


NMD 


' 


attaining  theraal  eqplHtxrlua  reported  by  Mesas,  Ma-far,  and  Pairbaa 


Coacluciona 

Uoine  a  tvo-particle  separation  of  the  -wave  function,  m  have 
aade  an  approxhaate  solution  for  the  ground  state  of  solid  helix©  three, 
y©  use  a  description  of  the  solid  in  a  sin&kj-particlo  representation 
vhieh  is  deteiinine&  ly  the  saaay-body  properties  of  the  solid.    V>q  ade- 
quately predict  the  density  behavior  of  the  cohesive  ©nerey  of  the  so 
in  the  body-centered  cubic  phase.    Boor  the  aelting  curve,  vo  predict  an 
antlferrojaogactic  Curie  teaaperature  of  the  order  of  0.01*  ll,  end  at  a 
sli$itly  higher  density,  a  Curie  temperature  of  about  0.003°  K,  in  agree- 
□cnt  vith  oxperiaoat. 


• 


ASOTDICEB 


.:,;i,  ■  m      ;      -.?  ;ion 


mas  Bamg/sgreaa  pammbe 

lie  use  the  following  vsll»4£nava  tfcMMU 

Iff   fix1)  -  f(:£x>  yx>  ax,  x2>  y2*  y^>  •  •  •  *%}>  pgi  n^,  rjgj 

*-->  .  .  .  Rjq*  r^,  .  .  .  i,IJ_1  g)    is  ft  f    q    v&rtc&lm, 

*nd 

u    "  "l>  *1*  vi»  ^c*  "  -1*0511 


i  i 

and  if  the    s      are  eecH  .tstftcne  of    a    TO?£*££es 


dvT3       i«l 

2n  our  case,  ve  have  a  spatial  \.wro  function    g(u  )  -  fi^u^  v^, 
v.,  IL.j  .  .  «  v«,  v„),    uhere    u,  v,    arid    v    arc  the  Cartesian  coordinates. 
She  nucsbor    n    is  therefore    311.  luvinoge  vave  function  1b  a  uev 

function    t(xX)  •  fds^  y^  »x,  l^j   .  *  ♦  y^,  Up  fjg  r^,   .   ♦   .  r      Wfa 


»   *    . 


K).    tie  noke  ouch    x      a  .function  of  the    « 


1 


X*"! 


ihe  raft  of  Pluvinage  and  |dtt  and  Dorovrits  displays  the  I 
placian  operator  in  cjfaerical  coordinates,    Mo  peffftim  a  sct*  -Priva- 

tion hero  fa  order  to  display  the  transformation  in  a  forei  iade.  of 

wmmimjiteiktu 

'in  our  notation,  the  upper  index  denotes  the  entire  class  of 


variables.    For  eseoplc,,    w    includes    ::, ,  y.- ,    and    s« ,    as  veil  as 


r,.  .     Shcrc  is  no  syahol    y*# 


ij.   .r..    at 


■ 


i  "WW »mimmb>mm<im+  Bvummt  »■■■»  m, -mm.mmmmmi'mmimmm  ■<  i  ini'W»» mini* ■»■ 

«•  /(ux -  xL^f  *  (v1  -  v2r •  (vx  -  v^r     (a-2) 

«tc» 
'.oaatruct    *      in  tl&e  v©y  so  that  the  operators  In    2    can  find  ap- 
propriate variables  in    f*     Efewiwrj  the  rules  of  stymies  ajgxly  to    g 
so  ve  cmst  retjuire  that    ?  a  g    cveryv&ezie. 

qpemtioiio  on    s    using  the  function    f ,    and  not  involving 

derivatives,  ps«ocaod  n»rely  trj  aprplyine  the  conditions  (A~2).     ncwcvc 

(i)  (^  )** 

if  derivatives  uvch  as    &g/^         and    3~c/*u  &se  required,  wa  aust 

CEsploy  (A~l),  then  apply  (A-£)  to  the  resultine  functions,  if  v©  viah  to 

use    t    no  an  analytical  tool. 

For  essogple,  vo  obtain    ^e/^v**    "ince    f  ■  e    everywhci- 

Sv      °    3v     " 


tban 


^L   »   |L   «     r   2L.  ^2si 


%  3-l^5     ~>v 


i:  "  &  *    *    §r  ^  (A^ 


}yi  <  >  k    —J*      "*i 


Shis  t$  the  rule  of  differentiation,  Eg..  UI-(6). 

So  foro  the  second  derivative,  vo  onploy  (A-3)  twice: 
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3v       5y  k    KI    5v. 
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■JjDO- 


£        1J  1        k 
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H 


^r 

g^,^^  »%•» 


JL      &r, 


it  ;.k 


5£r       5rll: 


k    J  "ik    ^ 


„      Iff      ^rlk 


*    ^  to 
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o 


^fj 


I      P*JI        **        I    «.,  ««»«i.«ii«|t»«»  (un   n  H«»l 


r±S**t 


4 


\^ 


v1  *""  '•"v*"*"    v^- 


S  y^f>mj«y  "■-■■       »  ,.  JTJ& 


— «4* 


o   ft 


V  111.   «.,!..  +•  '_,  .IP  II 


(A-U) 


Ih  %jh»t  follow,  wc  label  uoth  the    hl,  y  t  v,     end  f  y«  z± 

fcy    2T. ,  y.     oad    .L .    Ite  oo©&  only  to  roaaabe..-  thct  the  coc 

ftl 

1 
=«■     coraauto  vith  all  the  01  ,iol  oporetc. 


&* 


n 


o?  (a-4).   rutting  y*  •*%#  ona  <3c 


^% 


*  — -»    *  &  #    vo  obtain  the  taplaclc  ,      lace    r. 
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ulsj 


o 
5  r. 


LI 


*i  "  *<P*. 
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Putting  (A-5)  into  (A-fc), 


a 


£      {     2J^     «  a-* 


I      X     ^4^ 


r      '^  *  *  f*  T* 

•J 


(*I     "    ^ 


"I  J 


ij 


'ij 


here,  ve  have  separated  the  tern  involving    37—  <c~-    into  ito  cycootr.. 

cad  cuitSQyirwtrdc  parte  1 1  and    J.     Cu  perfowdfi  svctmtioa  over 

x    end    J*    wa  find  that  only  the  qpNVtxlc  port  rccma 

.  ;:loaover    ::}     aupprcaoilrjG  ?,     tnfl  that 

a 
-  jc.)*"  *  x\  (9    vo  obtain 


;fc    ^ 


«*' 


t 


-  u 


bU 


ij 


+    ESS  — «>-  -  ^i.   ,--£-  ♦  2  s  E  ";   — = « — x ■—  jcr-** 

Ife  introduce  tlie  tronofonaetion 

av  t «  •  [the  lnterpsrticla  HHynlf  umiit)) 

(the  position  of  the  cent*. 


TIM  1 


■  -  Q  « — 


1  a 


-  2  '7  )      , 


».JX: 


oo  that  8  (l^  -  Kj)  C:^  -  s^)  *  s,^  •  s^  . 

After  putting  these  SeiTiaitlooa  in  (7)>  ve  got 


-a. 


id     ij     "*ia  *4 

uhere    7,     and  <jrerte  on  fuactioao  of    Xj  y,  s,  o    and    t,    m 

where    «c~    acts  ooty  on  Amotions  of    r    •    Ihe  conBtmint    r. . 
i  i 

«  o. ,  •  i- ,     U  to  be  used  after  dli"  ..     Sberefcy, 

(A-G)  on    £ ,    the  »ec  I  the  scex»  go  t33ouc&  *w©  bo&  proceeded  throu^i 

usl^G  the  actual  function    g    wA  amy  apace  GKrivativeo.     8 

pslcseo  on  the  Gwctaotioos  indicate  that  oil  M  are  acto 

laeeher  of  the  nun. 


• 


•:  b 


JTXQH  .. 


She  iMinlifflUii  amplitude  to  Saflned  "br;  2q.  Il«(^6), 


h 


V(r)|x(r)«0 


*r2       r 


V(jt)  *  k  c 


r 


cr 


(5-3) 


viier©  \ie  drop  the  Gubocrl£rt« ,  1 4    £rca    r, .. 

ore  xmrticu2or2y  Interested  la  the  solution  In  t*ie  m-jQc 
r  fines  Just  inaide  cr  to  infinity.  An  i  -..on  off  X  in  powers 
o/r    vlll  converge  In  tule  nanco?  HI  (leEaonoti'ntod  l  ~l) 

in  detail,  \»  have 


If    q  »  cr/r,    then 


d2     ,    2    d  £   jt 

dr  ^         o      dg. 


g 

Ml 

I      / 


X  «  0 


->:- 


oy  a  otrai^tforunrd  trtaiofonaatlon  o£  doxlvatives.     Shen  (k)  baa  the 
fora 


>   a' 


a  or  kth 


a     — «— ?j 
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.  qu    X 
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2hc  vnrienle  q  xrl'U.  not  vaaioh  la  the  range  of  interest  so  v© 

B 

divide  lay  q  .  tto  eloo  define  B  ■  iL- — y5^  ,  ooa  factor  eoateinlag 

h^ 

k 
ell  t&e  Bhyeical  peraraotero.    Por    if*,  B  »  l6.60j  for  I3e  ,  3  °  .  

Vte  nuot  solve 

g    ■    B[qG  -  q2)X.  (B-7) 

*fe  A  solution: 

CI  q  *  C2  *2  *  c3  S?  *  •    *    • 

(B-Q) 

Putting  (3-0)  Into  (B-7),  we  obtain  conditions  on  the  eonov 

q**1:  0*30       •  3  C  ^  0 ,  «  0 

0  *  Cj 
q°:  ^-IC^-IC^, 

i:     HCrDC-rDC0  , 
q3:     5*C5  *  BC^  -  BCX  , 

and  oo  forth. 

y©  expect  the  solution  to  approach  *  constant  val;  tpwin 

of    q    at  lar©e  !/<!•    tieaao,  oil  the  Lest*  with  no&&.  \&mrt&* 


x*  afcArtt 


• 


oust  vanish  (in  other  vov&a,  ve  hope  to  find  o  function  cotisfyine 

1  2 

vitfc    x  *•  cQ  +  C,  ti    *  C-,  (    ♦  .  .   .  j  uii&out  negative  powers  of  <i). 

fet&ng  oil  the  C's  vlth  negative  Indices  equal  to  zero,  ve 
cfotr. 

c0  *  C0  C6B° 

C*  «*  C\  cl  **  o 

1        1  7 

0>  -0  CQ  •"C0B?/(6'7  *  ***3) 

c3  «  o  «  cxir/(--3  *  5«») 

C^  «  -  B  CqA-3  C^  «  cyi/CUK;) 

C5  «  -  B  Cj/5«%  Cu  -  CjS/di-lO) 

u12  32-n    v8 


ciuM  ]Ar<V^ 


'13 

<W   B&K   (ci--cib5 
%*  "  5P2g   ^cifc.  "  *W 

Ve  shall  demonstrate  "below  that  the  solution  conylng  the 
cient  C±    is  non-xtoysloai,  Hence  ye  do  not  display  aany  terns  of  that 
solution. 
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We  use  the  valus  of    B    for  tier  (B  *  36.60),  aa<l  obtain  mow-. 
values  of  t&e  coefficients  ao  follows: 

C^  •  -  1.391 

:8 


CA  «P      &.1?  x  ID" 


C^  *     1.355  2  Xf X 

C12  *  "  5*25  "  X 
c^  •  -  1-^30  M  2D-1 

C16  "      3'6:5  *  *°~3 
3.05  M  M 

c^  *     Ml  x  10"3 

Vto  can  write  the  •olwbioao  ae  foliovc; 

x»c0?0(ct)  *cxar1(q) 

PQ  »  1  -  3/<4*3)  <t  *  Ba/(8*7-4.3)  q5  ♦  B/(J0-y)  q"^  ♦  ■  .  .  . 
Fx  »  1  -  B/(5**0  q*  *  ^/(9-Q#5'«)  C  -  3/{ll-20)  q20  +  .   .    . 

At  large    **,    the  phsyBical  ci  .racopoaAft  to  two  helivc 

atone  floating  in  apace,  with  ao  tendeney  for  re  -  action.    ttoj  there- 

fore, eaptot  the  total  flow  in  or  out  of  a  ejsbere  at  a  large  distance 
be  veaiofcingly  ooaii. 

To  a  low  consistent  order  in    q, 

X(r)    -*  cA  -  c,  a    ♦  c,  q. 
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-    J* 


2btalFlus«    Ha    J  dG  [co  - -^    ♦^f  ]  Um    "[Cq-^ 


C,  ok 


r  -•«  i? 


_  %  ih  Crt  C.  c? 

*  .  .  .]  « -■  ■■■■-■  .  (I 


HQw,  if    C0  •  0,    then    X(»)  -»  0,    in  violation  of  the  condition 

x   not  offset  the  particles  at  laxee  distances**    If  CQ  ^  0,    and  if  the 

flux  is  to  vanish*  ye  have  the  conditioxi 


^If  us  have  a  bounded  system  end  the    $4 (r,  5    venis!:. 
finite  distance,  then  the  aeauaajfcion  that    Hl.<ri.  ~*  *>)  *  1    loses  its 

aeening,  since  the  toovsadaxy  conditions  on  the    $.     will  prevent  large 
values  of    r^ .    frees  occurring.    However,  for  ays-tees©  vltfa  Atom  ■ 

such  that    V  >  >    o,    we  can  aesuae  that  m  have  m  unbounded  oyo 

in  setting  *$>  the  conditions  on    X#    Shis  condition  io  :"iot  restrictive 
for  the  ras^-bc&y  cese. 


In  the  theory  of  Chapter  IV,  uc  have  need  of  a  non-local  average 

2 
of  K,  $    the  square  of  the    elation  function,  mid  of  JL*  "7  X4  4  , 

a  correlation  eneitjy  density.  In  this  appendix  ve  derive  a  (rencr-til  re- 

oult  for  the  integral 

H  5*1g«  *  ***  *o»  the  coordinate*  used  in  evaluating  (C«l). 
Particle  A,  baa  a  diaplaceraent  q  fraa  a  lattice  point  J.  3»  particle 
A-  ia  deecrihed  ly  a  probability'  deaeity  proportional  to  e   ^3  .  2he 
function  ^r.,)  i8  **■**  function  xfcich  io  deteroiaed  tiy  r,.,  the  die* 
tance  of  A,  frca  the  field  point  1.  It  could  he,  for  exenple*  a 
potential  energy  or  a  correlation  function.  Ub  eeaune  that  F  is  a 
function  with  a  cut  off  at  r^; 

or  ?  »  ^i^  8(ri^  •  ~G*" 

In  the  second  e;cpreeeion  ia  (C-2),  wo  have  cs^Ooyed  the  step  lea  to 

provide  the  cut  o: 

She  function    *(*<«)    ll  Civen  by  a  aeriea  eaqpaaaioa  in    tj/r, 
vhich.  does  not  converge  for    r  <  e  .    She  function    XO&j)    Ml  fcttowa 
decrease  continuously  fron    r  *  **    and    >'.{c)    is  very  eoall,  so  the  i 
ditioa  (C-2)  ia  a  {pod  description  of    X    for    r  <  8*3  A,    and  eliminates 
the  problem  of  aco         Wi'ltfMii  of    Ha/r)    in  that  region* 
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Figure  12.     Coordinates  for  Long  Range  Averaging  Integral 
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Qtm  wo  vill  find  that    F(ra  J     la  siow^  varying  naev    r.  .  *  u^, 
vo  vill  eagtond    *fo  .}    about  the  point    j.     2»  resulting  eapraesiou  wi.ll 
tic  a  function  of    u^  »>.    the  distance  froa  the  lattice  pc&at    3    to  t!& 
field  f  !•     Using  the  susd-lia;*/  vxiraablec    a    end    p,    to  tbe  order 

shewn,  we  have 


F(r)  »  F(u  -  a  -  ^.r.t.ft,„  } 


*  F(u)  -  (a  ♦  e~^~^}?!  *  -'  ...     . 


Ito  integrate  the  first  tuo  terns  of  the  tagftMtafc    Ftmt,  v&  beve 
I,  «  F(u)/  e"6^    e(r  •  rn)r2  dr  d(coo  e)  2rr  ~: 

•  »  F(u)  /  .-  ^  ♦  rS  -  2^  CO.  9)^.  .  rQ)ra  te  d  ,,.  ,  jL  . 


*•*-»» 


Bo,  tie  est 


Wo 


U  D 


-d2(r0  -  u)        -Q2(r0  ♦  u) 

o  -  0 


♦  ?<u)  f  tl  -  *  [5(r0  -  u)]  *  I  -  $  [5<r0  ♦  u}]}. 


Ucw,  wo  vill  find  that    d    >  1,  rn  »  2.5,    so  that    rA  *  u   e«a  be 


■0 


0 


eidered  to  be  very  !&.&*#  end  \m 
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I  [1  ^<f>  [6(u-  p0)]]         C- 


«&2(r0  «  u)2  - 


?8u  J 


Me  ne:;t  com.  o  torn  contr  F'(u).    Sioao  i*»  ac&toetioa  tcn» 

with    t&(r0  *  u)]    la  crvuluotipc  the  lec&inc  torn,  we  ewe  Justified 
doing  so  k  .  oct  here*  sod  ctf&li&tc 


<7      Un»j 


0 


•  F'(u) 


a  i  ° 


\ 


!3mh  » 


^  •  •  r 


-Jj    tX  ♦'♦Wu  -  rJM  ♦  ~ll^-  {3  ♦ 


*n) 


X    c 

Hi  have  the  rer- 


/  •"**   F(r)  0(r  -  r0)dv 


/  •"*"« 


■ 

28     uj 


I 


-  P'(u)(u  -  r0)(3  *  •  vQY 


(«  -  rQ)3. 


|tt  *  Ij 


♦  0  (F'Xu)), 


-112. 


IgMfcUai  (C-3)  cua  be  use&  Ic  tbct  Will  na- 

tion* used  da  evaluating  ttse  i  y  •«©  vlll  vase  tbe 


Awaanax  d 


—  m>im.i^ii     !■■■■■ ■»»mmi»'«*i^M»iii'.»..Pii».ii     »i  .,  ».|,pi  imi     ■■■inw»wp>u.WWHnii»Min-.iwiu     >,««■■       W 
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^1 

<h." 
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-£.     **      ** 

1 
d 

^1  *  ^  '  "*ere 


D-X) 


ftpott  Eq.  W-(I7).  t,    wbm 


f  -  Hi  ♦  |<oU  -  rn))3 


D-   • 


+  0 


2 


"7  -5w(u  -  rQ)' 


For  valuer;  0*    u    vfel  ronalder,  the  second  terra  io  alueys  much  eoaller 

than  the  first,  so  that  the  denominator  of  f '/-    cn»  be  token  to  be  the 
leedine  tern 

*fe  define  a),  f    «  |fl  *  ^(ftfu  -  rj)].     Excluding    r 


of  order    -9-  ,    ve  imve 
6  u 


f »    ~ 


■ 


f\ 


"  g"       ™»       1  j  ■  1  'jy" •        *»       ""y       *       ■  "■■a*"      • 


'Ihe  first  of  these  la  (l>l)  gives  ga-j  the  second  contribute*  to  g: .,. 
Una  correction  tern,  A.  ie  of  higher  order  la  1/6  and  X',  end  accord- 
:;ly  neglected.  First,  ve  use  the  series  expansion  representation  _  r 
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«    ^  Ji  *.      2£  ^3  22 
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3-^TT    HT^'M   ~7 


fcy  difftere  (j>^)>  and  tken  d:  v$»  a&ct  : 

[X/u  ♦  1/2  d/du]  f'/f,     getting 
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Ac  a  (&*$)*  w  an  c.  fci  Yi&cwb  rtpAtsoeirtcitioa 

far  -  X  ~  (r^r)  *to(«9<y  *  r0)). 
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t*e  obtain  the  cor  »^    "by  ^itt|  out    t*     to  second 

IT!     Fxm  the  HMBlXt  of  Unjf   .  Ml  oall     f  *«j)l  then 


■flgfii  . 

tx  «  (const)  /  e*      *    ^(^j)^* 
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Exclu&ing  torsas  ftp 

have  for  the  toroa  in    *     and.  tlio  lacaing  te«ao  in    X'X 


(u  -  vQ)    end  p    ve 
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(6a(u-*0)2)  u-r. 

c  *  ooaXJUtr  tc 
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Xa  arriving  *fc  *&U  I  fern*  ve  fiad  that  all  the  eontrilju- 

tlcos  in    o    to  the  load  f '    have  eaneellel  with  teras  fr 

the  oceoed  tc  !•  correction  tern  la    0    arises  -  -c  owl 

off  la    rf%     Shore  li  ftlaO  a  collection  fcers  arising  r*oa  the  slope 
in  the  region  'che  Ctettsoiiia    ejqpC-  &    a.  )     lo  iarec,  but  thio 

/o4**  compared  to    Jr,    Wl  to  have  aiocerded  a 


tera  is  of  order 


tern  or  this  ordei*  :1a  arriving  -7)*    Wt  will  use  the  correction  beta 

.'},  but  cent,  ortolan .  -•  las  of  the  order  of  the 

change  in  reeulte  brought  abo-..  =>7). 


■  • '.::::", '..  ■• 


'•^cuiatioii  or 
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I&  Illustrate  the  tsettood  used  in  constructing  Pig.  6  ty  evaluat- 
ing the  points  at    a  *  3- !     - 

Met  first  evaluate  the  sua    8  - y*-      using  the  series  e::penBiQn 

fora  of  W-(37)  plotten  in  Hg.  }, 
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J  k  " 
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W*  uo®  the  Ginc^fi  'particle  fey  *acttrix  to  -jqpstHient  oil  the 

particles  i  p    and    k: 


- 


<  -^  A^<r*[s15<p»,,*i««s<k,^H 


a 


fcJV^k'n  «t 


(*4) 


La  (?~4)?  t&e  product  over    t    U  r^a-ditt&otsal  eo  tfaet    IT^    raay  operate. 

tfe  new  present  an  argtmnt  tlmt  tbe  miaerotc.  P-4)  vanishes 

to  lowest  order  ta  the  oncrey,  in  a  Braveis  lattice,     lot  us  oca 
etons  1,  2,  and  3  arranged  on  lattice  oitc  (  and  3  in  on©  o:*  the 

lines  t&j:  1  Iktsvais  lattice,     vlo  have  for  their  caatr 

to    ZX\     fron  tho  •     .a    IIe 


■to 


"Com*/  a^c  (  '■   ..  M  eU 


♦  Const  /  &;  -  JU> 

Hero    1^(^)1  ■  x(rjg)  x,fri2^plP'    ^ 


3  -5) 


P3f*l>J 


4>(V4>  ' 


la  \F-5),  *t  have  put  tlie    1    M&i  alone  tho  1:   1  g  tbe 

lattice  cites,  and  aoouac  cylindrical  syix*.  to  that    y.  -  |L    and 

a.  -  a,,    average  10  zero. 


■u  ■»". 


I 


• 


.  tie  tfM 

■ 

the  oe»e  pe: 

3 

•   )  fey  vrritlafc  thr 


a 


*  /  d*J-  a  ♦  ad  - 


] 


)  - 

Ixx    f*&),  «•  hove  .  time  OB  the  WMMp 

- 

9ft 

:  •  ■      • 

Oft. 


•    I  .u 
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A 

• 

.•at  relationship  stems  STroa  the  fast  that  where    g{l,2)    has  an 
appreciable  value,  near    x^  *  °#    «(1*2)    is  aoaci&et  taore  gently  varying 
than    e"    *  .     T3m  second  is  cotained  tjy  inepa  .  Pig.     i 


Slierefore,  ve  neve 


H(l) 


£    - 


h    , 


I 


1  IT    -1   A 


4s( 


f  fcfe*6!  '12  ■ 


Iherc  are  h  such  pairs.     Hone©    U*    2$    1*6° 

In  a.  Gijailai-  fashion,  ve  can  eetiiaate  tlie  second  te^n  of  -F-l). 
Because  the  atcos  aix?  localised  near  lattice  sites  in  en  ordered  lattice, 
but  the  suscatloB  over    i:    yields  only  the  single  tena  where    sV^  *  -  e... 
Putting    1  *  1,  J  «  •  %    end  using  t&e  sane  ftpprosdaatioaiD  as  in 

(F~7),  ve  01 


^2) 


ta    «. 


J  22 


ir    l  r  '^12  *  *i3 

5   I  '    r^r^     »      ^  X12  A13 


x      n 


if     ^W^ 


•liV-tHi 


3t 


of;,  di 


.  : -: 


VH 


-iaa« 


m  ■■■■ 


^  (U)  (.01)  * 


'iftei-e  cacti  four  such  poire  per  atom,  oo  that    IL'  v  «*  3-2°  K» 

She  eetimtinc  procedure  weed  here  deaoastrateo  only  that    H'     le 
not  tbe  tadaeat  tern  in  the  energy,     3ae  assuaptian  tbat    H*     io  onail 
end  has  little  effect  on  the  uave  function  isust  reraaia  oa  aamsaptloa  of 
tbio  theory. 
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la  this  appendix  we  obtain  the  factor  l/3  Introduced  Into  the 
energy  in  Eg..  I7*(^7)  by  c  £*>rc  detailed  otody.     In  Eq.  Xf~(2T)  the 

oatrix  oleiaeat    (IT  1*7  >  ft  )    ©ensures  a  "klnc 

t  .  X       J       fc*\3 

of  particle  l,  ' . ,    the  part  of  the  wave  function  of  particle  1 

■which  is  determined  by  correlations  with  particle  J.     H»  sub  of  all 

such  cant  the  sua  total  of  the  single  particle  energy  useult- 

|  from  the  caafiaii-  NMh     In  imtliecsatical  sysiboia,    rA  Jf. 

hS  n- 

I*t  us  eet  up  a  coordinate  frcwe  with  x.  parallel  to  .L  ,  so 

- —  i 

thcx     'f  »  -  h*"/sn  d  /)x/"    -     ^et  us  further  easuao  that  Mr  syi> 

y  io  worth  that  there  are  two  jaore  atcos,  say  k  and  t,  at  the 

distance  11  s    but  on  the  y.  and  s..  a-«s.  She  kinetic  energy  osao- 

;ed  with  correlations  of  1  with  these  three  atoras  is  then 


• 


gfi   , 


<*>i^i = r  fc  [         -j  jj  -" 


f_  *•     -^-jr  n  f 


R 


flf   -*-"??  f& 


't'  v  2    t   # 


]] 


(0-1) 


con  write  the  sesae  c  r  two  other  sets  of 

x  -*  y,  and  x  ~»  s,  add  the  three  results,  and  divide  by  3. 


■  ■ 


V/i,< 


u 

at 


*.* 
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s2  . 


A  * 


(I 
1  *  "2  *  - 


V 


t    aft      n 


But  ia  our  iHlglllipHrtiliitl  repsroaentation,  all  the    f 'c    for 
tiw  «t  tte  am  dl«faK»    :%    M»«»  Jy  aytatalc  about  tte 

gin,  lat-  -to  1.     *!»  I  pts  an  the  derivatives  ore  then  oujx 

fluouo.    Hi  can,  therefore,  Witt*  £br  (0*2), 


<*> 


B    3 


i  ^ 


^    of  J 


rffe 


+   J'.:  7f 


/ 


a 


- 1  ^  ^T7^  * !?  *W 


(G-3) 


Writing  (0-3)  for  ell  the  atone  in  all  the  rfcqp,  one 

ve  bave  0  q-,fc'  6' 

,      h*    1    Jjk     - 


i.      It    1    _fi      - 

1— &-J.. _ 


dv 


CG-10 


/  o  dv 


1 


> 


W&  can  understand  the  factor  x/$  in  another  uay,  by  applying 
cXacaiccl  laschenics  and  9         ftMfcUfetng  the  dyaesical  variables  tttt 

the  corresponding  qaaatm^aocZioaical  operators. 

Xn  the  Oiagran,  ve  have  q., , 
Mwurtae  tte  aupiK«Ht    uj 

the  particle  1  from  the  lattice 
^  cites    3,    by  the  relationship 

u.  •  %  ~  Km  *,    or  eauivnleatly, 

Uj8  -  %  *  «^  14  cos  9y 

(C-5) 


We  vioh  to  sao&ouro  the       cnercy  of  the  particle  l  by  ob- 
serving tt*c  lengths  ot  all  the  u.  and  tehing  a  suitable  ovcrccjo. 
us  aaauiae  that  v&  cannot  sioanure  our  absolute  pooition  in  space >  but 

nereiy  can  aeacuro  the  scalar  quantities  u4  «  (She  probleu  ve  vAsh  to 

3 

study  is  the  inverse  one,  vherc  -^c  noasure  quantities  to  terns  of    q, 

and  deduce  the  corrcspondii:.  in    u..) 

\te  pr Q  ;oed  by  to^inc  the  tine  deri-vat:  ve  of  (G-5)  and  savoring: 

2qq  "  2q  V  COB  6,  -  os  ©  *  2  u,  tt, 

P  cos  e,  «•  c  ^oaX)2  -  «rf  1  coa  e 

cos  ©,  *  8  q  q  IT  cos  9,  cos  9,  •  u,^  u '*   . 

£uuad.&g  owr    ,3    and  averaginc    u0  v,    and    q    over  tisae,  and 
ofooerviae  tfcift    u,    i»  nearly  constant  and  equal  t  we  {jet 
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^  <T  <f  ♦  Oj  q~  y  I  •  q-  IT  n^  \coo  0^/     -    8  H^ 


J 


Divide  by    H,   : 


\  12 


For  aaoll    q,    or  lew  angular  ncaenUBj  ve  e©^ 

a,  '<t  h    *    *  *«2 


2 

9 


2k         it&d  energy  is  fchrae  tinea  the  overage  of  the  "kinetic  eaergiea" 
relative  to  the        [  I  polato. 

la  our  car>o,  vo  asii  for  the  apparoat  ■  hMllil  energy  of  relative 
notion  to  a  eet  of  pointo    J,    in  temc  oi*  the  kinetic  energy  of  a  particle 
.ag  in  0  path  deocribed  by    q, 

*         U.    »    JL   J  .  .  Cc-6) 

3  d  Q4    J 

2y  the  oarreepoodence  principle,  the  ritfit  lUi  Of    ■ -6)  repreooats  the 
-ective  fcine-  #*Wf*  •*  Eq»  XV-(26),  and  the  left  side  repreocnte 

the  nodif  energy  of  IV-(27)«    Tfc  are,  therefore,  led  to  the 

factor  l/3 . 


la  thie  appcndi::  «o  attend  to  the  detailed  problea  of  narnaliaiix* 
the  ground  atatc  vave  function. 

In  Chester  2V,  vo  found  that  the  eKpreesion 


„ 


r(D  -as  lid)!'-   8  /  rttO  **(raM{| 


could  be  HI  in  the  fora     r(l)  "  £  r*(3-)j    tfsere 


ri(l)«e^^   e    °    *    . 


Here    o^*  ■  -  --57  Ai  I!  /r('&)  nL{*te)^L'    &  •■■HHF  •which  ve  can  obtain 

k 

by  a  In.  sua,  and 

<■» 

ri-\ 

.^teet^ate  *  dete^ne  .to  effect  on    ^(l),  and 

this  neoeoeity  leado  us  lato  the  Mi         viion  of    A,    and  the  Hw III  11*1 
1  of  the     fyy 

;e    <j>. (*4)    ore  a  oct  one  vhich  are  elQdnetatoa  of  our 

si 

q  order  ISmLitonian.    3facoe  otateo  ore  occupied  only  because  ^©  cruet 
cjfy  tlie  Pnuli  principle. 

I  the  variable    I.    doeo  not  appear  In    IL,    tfoe  labels    i 
on    ^<(-m)    •*•  cood  qpflatua  mxdbare>  m  a  sjoy  •M0t)gf  parallel  to 


•1S7* 


i  M 

■ 
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tbe  treatoent  of  the  aboil  aatlel,  we  call  upon  tJicsa  fee  doeoribe  ti» 
atatea  pesuitted  by  the  I'5©!!!!  principle.    Seen  of  these  states  ',*?  w.,11 
assuoe  to  bo  occupied  'by  two,  ono  .  ■».     (v7e  label  tiic  9HMM 

■by    i    and  cot  two  particles  gfl  -ur  codified  scbeoo.  rather 

t&an  by    ,1  »  i  *  l/2,    fear  one  pr.         •  per  state  ia  the  visual  echcsae.) 
Slieu,  since  the  particles  arc  "..opt  fa  CtM  ty  gf  l«&t£o*  c.too, 

each  VMM  bo  or  ft  wavelength  attti  to  apaa  two  lattice  citce. 

la  Chapter  IV,  ye  wrote  a  pttttftl  set  of  waves  for  a  oir^le  cu. 
lattice 


COG    EJS  +  ■  coo 

yk        la^a 


2UC  *  g 


f 


£        vole  chosen  for  each  wave  to  pub  ao  nodes  at  2c 
ni 
and  for    n    or    i/a  odd,  to  place  laasclm  at  lattice  points,  or  for    B    or 

l/a  oven,  naxioa  oidwoy  between  lattice  point;;. 

I  will  00V  evaluate    A  I  particular  case,  where  the  lattice 

point    XL     ic  at  the  center  of  a  coiid    (Hv  a)    by    (U    a)    by    (U    a). 

•~  •*  «* 

In  this  cose, 

(H-3) 


1*    ■ 

coo    -r-~ 
n  a 

°       W*Xtf,t 

u 

a   *    j    ■ 

»*  «  x#  St  J»  *  .  *  -  (n)^3* 


evaluate    A    by  writing 


.-.,:•  v  u  MO 


/  rxU)  4l  -  a  /  s  l$/a)r  8  /  ffcOO  x 


-w 


*- 


;n-i») 


5W  fa. 

0  % 
We  hove  leaned  that    *  is  fast-v&ryi13®  cenpawd  to  the  lat/. 

apocrine.    Benee,  ail  the    cos      tenao  can  bo  taken  "to  bo  unity,  end  we 


or 


/  r-jU)*^  *  i  -  a 


i,£ 


A  * 


SJl  c^ 


t*e  now  use  IV- (31*  5  to  eotiaate    a 

-d 


Ua   a  t  Mu  (*.)(*  -  o 

^     i     a    ■*■ 


lite  will  eatiante    a    for  the  aeso  leading  to  the  wave. 
Using  the  lowest  terras  in  the  cc-gpoaoion  of    coa  (x),    wo  cot 


**r% 


i 


s  c;      i^c 


"      u^t,y,s 


1  - 


wvv. 


2 
u    it 

" 
2ou 


W*  cannot  evaluate  the  auaa  in  (H-6)  unless  we  know  none  about 
the  relative  va  «f  the    C.^.    Ife  now  show  '  f  we  use  the  assusp- 

tloaa  pre-  )  takes  on  a  parti  ;;  siia&le  form  because 

the    C4^.    are  all  of  the  saoe  raagaitudc. 


-::.;P- 


fOXBJ- 


l*e  write  the    9    body  denaity  raatris  in  a  slightly  different 


r(B)  »  i  pc  I^CDI2  tV^)  ^  vu  •  •  -I  N^ttf 


p 


2 


1<W  XS3  'S*  •  '  "]  '*3(3)r  fV'r3l:  3fe  •  •  J 


,o 


•  •  •  •  'VH)        a^»  •  *  '  ***>£•  -7) 

la  (U-7h  vo  eooin  neglect  overlap,  tfe,  therefore,  have 

Nt«Mi  tx3i  •  •  •  V  '  •  •  'VB)I2 

Because  MMta  state  is  occupied  to  the  mm  ei-rtent  as  any  other  state, 
ia  the  woy  we  write  (IJ-7)  and  (E-3),  all  the  pemutationa  cive  the  ssras 
value.     In  (B-8)  vo  have  written  n  particular  one,  and  aultiplled  by  the 
nuciber  of  poraiitatioao, 

W&  Jaarnod  in  Chapter  IV  that  the  correJatlone  functions    X 
con  be  treated  to  a  good  approjdiaation  ©a  havine  one  of  the  taaabero  of 
the  pair  fastened  down  to  a  lattice  Sit*.    'Hum  the  product  tX^X^JLjL 
.   .    .  JC-J     li  the  imam  root  of  the  :„■         ■  product    H  Xr  (tt«v)*    ufe 
we  evaluated  la  Chapter  2V  as 


i 


Wo  con  \*xlte    ?      la  a  fom  vhicfa  i»  general,  for  the  oat  lid,  p 

via<*i  that   a^  >  >  4>  : 

at 
I 

s         a. 

Mt  only  one  aoaber  of  the  sin  has  cry  appreciable  value,  becaxise  of 

our  mempblm  of  no  overlap,  ye  can  Mi 

ft 

r-  o         \fc      - I   r 

'.    »SE  Dr(a)l  *        e 

Me,  therororo,  can  vritt  (H-G)  in  the  following  fona 

Iter©,  we  have  Mktftt  advantage  of  the  .ily  the  eiaelG  nec&er  of 

the  sua  over    J,    -uhere    TL  «  «L   ,    giveo  any  contribution.    Each  of  the 
etateo  equally  populated  with  two  porta c Ice,  hence    fc  «  1,  2  .   .   . 

and  all  the    I    integralo  are  equal.     Because  of  the  fast  varying 
Oauaaian  function,  MOB  Of  the  cjGinots  in    $«..!)    io  practically  unity 
in  thie  integral.     '-Therefore,  we  have 

o     -*       ■*■ 

ft 

Define    8  51      *  (3(a) 

V  •    V  C^    «*)    /If  " 

0 


*JC^»*tt7      '". 


or 


S03c 


Cance  ail  the    (  ore  equal,  Ml  flftta  t#    (&Q3  o/[(2::)3'2  I  p(a)l)1^ 


n    cQ.    Ike  numerical  value  Otf    GL    vill  not  ue  useful  to  ua.     Ifcwever, 

the  fact  that  the    C. ,,     ere  all  equal,  vlthln  our  QoaxaagptloiiB,  is  of 

•  - 


value. 


Ua  1—1 1 tH  (H-6),  mtoat    Cq    for  all  tfet    C4#,*t 


2 


i       n 
c0m    " 


i.    •»   *""X      ""*"X 

tin     \ 


g 


(H-UO) 


an  the  2ee&in(s  cams  la  the  product.     (Stje  o:!panslou  Hi  use 

(H«&>)  la  eiipeeted  to  lie  valid,  since    i      *     will  provide  a  cut  off 

In  the  region  o  ca-iveryenco. ) 

fit  the  sua  of  (H-10),  thexx*  are    2(11)"' J    tenas  with    ft.  »  1, 

:i)  ,:>    tenaa  A,  "  3>  -ate.     «bia  cooes  about  hecauao  the  quantities 

l/*W  *A**#  -/%    oa^o  propoiTtioual  to  the  ccqponents  of  a  vector  in  the 

reciprocal  la  .  of  our  aeauoea  Gix^ple  cubic  lattice,     Sher©  are 

ja/a 
.    •  (t?)    "  x>intu  .e  plane    (L   *  (coast)    and  an 

equal  xns&e?  in  'die  reflected  plane. 

wa  can  there      i    rewrite  (3*10) 

-cTaf 


X    «      Sst 

if    I  "*        "  O 


a 


1  ''       » * 


*♦  I    w 


I 


*/*(!♦    I — I)      »($♦...  ] 


-.3::. 
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Bat    ~  »  4r  ♦  -4  •  .  -  .  **   f  *  4?  *  •&  ♦  *  . 

r    2** 


2* 


•111  "y^ga? 


.2 


a   » 


for  largo    V«     ISeooe 


4  JB 


I 


<f.         1 


.i 


a 


|r    (:r  ♦  y~  +  I 


/j 


■  £ 


o~^*7 


or 


>— 

■u 


fn-n) 


In  Bq.  IV-(55)#  we  stated  that  the  contribution  to  the  enoiigy 


4Ti 


P  /  ?*i««V*i(1)*        * 
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-    1(0     ..   ._     -•.--»»: 


.,  c 


(B-; 


and  la  of  orcter    m^***    We  now  ahoy  this  explicitly.     Using  the  roprc- 
eorctatlon  (B-2),  wo  Iicivc 


<•>• 


1       ,   co2*2 


IT 
a*" 


l  ij-  \  a,  /  J 


1 


— -^.  ♦ 


I 

A. 


I  0O°  I  rf  cos"   rf    e  -'    • 

in,a  /  i  a.  a  x 


wber©  wo  hove  act  tin  integral  in  tfca  tattQMtai  oqwl  to    A^. 


>  & 


As  '»?©  did  la  (B«#)j  VB  am  tafco  all  the  eoGiines  to  bo  unity. 
Ms  then  set: 


9      . .«:  -     feia 

a  a"  S    (l) 


A^nln,  the  eua  4S>,    &a»    21    tenaa.     Using    l/l  ♦  1/3*  +  1/  .   „  , 


i& 


*  st^/D,    we  hew 


-1    £   r^^/3   3Ba 


9 


oaoqo: 

Single  Particle  Orbital 
B.P.  8.  atraar<Jee  and  H.  Prisa&off ,  ftiys.  ,tev.  |]gj  .9&>) 

r4  33*0  position  coordinates  of  particle    i,    not  Including 

*  v. 

spin    £,   (3L     includes  spin) 

■She  caardinateo  of  particle    tj    including  opia    £, 
<3oe©  not  Include  spin) 

Decrees  Kelvin 

Anfiotroa  unit 

S  P  later' <aetoroaaont  of  \4mt  follow*    2»  penautation 

p 

operator  for  Stem!  partlclsa 

I  |L  Sbxxas  a  oyooetric  percttitation  of  vtaet  follows. 

P 

£  2*4  Paaautoa  o:i!y  the  coordinate  indices    1*    in  tbe  anti- 

P    x 

oyiaae  tease. 

o(x)  Qtroc  B  fuacti 

o  3bo    vnva  nudber,!;  of,  the  single  particle  represent© fcion 

I)  «  const  o  ,    q,-^  »  distance  fro»  lattice  site  1. 

rA  X|x    x     ...  x  )        -iaelty  aatrix  of  order  (p)  'nee  Chapter 

IV.) 

x  1±  aity  natrix  of  order  1.  (See  Chapter  IV.) 

a^.        S5sc  displacement  between  particle  I  and  particle  ,}. 

variable  has  a  derivative  ^z,  but  no  derivative  c)/Y 
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i*  g  .-aerpcuticlo  parameter,  union  is  an  independent  vari- 

able  in  the  .'^muiation  of  t&e  Scftrodlngjer 

equation.     It  boa  tbe  constraint 


Sd 


V 


mmummmm 


*u  '  /  *  (yi  -  *i>  +  (si  *  v  °"8^  "  0u  ' 

t&ich  jawot  be  ixsgposod  %fcen  evaluating  wafcria:  elenenwQ  c^ 
talning    r. ..    Shore  is  a  derivative  of    r4  ,,  3/V  ,  ..    bv 

no  derivative    d/!3:c. 

t..  33ie  position  of  tlie  ,:outer~of«iaaso  of  peptic  I  and    ,}. 

Sbere  lo  a  derivative    %/cte,  but  no  derivative    I)/V. .. 

X1.(pff .)  In©  correlation  function,  defined  by  fornula  IH-{17).     It 

hoc  a  series  representation    FJo/r.. ,)  &(***  -  pq)>    ^c 

0  ft§  the  ©tap  function,  and  a  Yukava-type  repreac  .  ■       n, 

-  r(r  -  p0) 

1  -  r_/r    e  . 


2be  tUstaao*  from  the  field  point  k  to  the  lattice  point  j 


^(XJ       xe  oing3e~particle  wve  function,  a  function  of  space  and 


"4 


..•  )  2be  ffUglfll  |)W1I  ll<  -wave  function,  a  function  of  space  on3y. 

•iff     2:      +2 
0U)  liio  error  integral    #2)  «  (2/st  '*")[    e        dt. 

F.  Hf,(«L,),    lunere    t    is  the  portion  of  an  atata  near 

lattice  9 

F0(jc)  .00  ecproocion  for    X,    the  correlation  function. 

.)  :-tsry  function  of    (r,  ,). 

km* 


Aj  k.^|  -  \\Ji     (^FiG.  h.) 


I 


«n 


i    f    1    4 


Q»  distance  feara  the  lattice  point  Oenotinc  the  ainglc- 
portlol©  ijtflWiBffli  pair/:;  to  the  alsell  of  atona  at  distance 
&i    ©.  function  :c  esecewtsy*       "cc  Fi£.  4.) 


Wro*  I2odine  Clin* 
Mnkhm 

YQxfeshirc  0-235? 


